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CHAPTER 1 


INTRODUCTION 


With the recent trend toward digital processing of signals 
and the advances in digital circuit technology, various digital designs 
of phase-locked loops have emerged. At recent conferences Cl I C 2 I 
entire technical sessions have been devoted exclusively to digital 
phase-locked loops. 

The development of digital phase locked-loops has proceeded 
similarly to the analog phase-locked loop development in that two basic 
problems have been considered: tracking a carrier signal (or synchro- 
nizing bit streams) and demodulating FM signals. Typical problems as- 
sociated with carrier tracking are the mean time to gain lock, steady - 
state phase error probability density, mean time to lose lock, and the 
location of threshold. Analyses concerning FM demodulation concern 
the loop performance based of signal-to noise calculations and the lo- 
cation of threshold with a modulated input. 

One method used to track a carrier involves varying the sampling 
time [3 3 in an attempt to lock on to the zero crossings. Another 
[4 3 uses in-phase sampling and quadrature-phase sampling in order to 
determine phase error. In addition, these sampling rates are variable. 
Other techniques [5 ] , C6 ] , which are not completely digital, require 
mid-phase and in-phase integrations, followed by a multiplication in 
order to obtain phase error information. 

One of tfc e earliest schemes for constructing a digital phase-locked 
loop for FM demodulation [7 ] employed a voltage c ontrolled oscillator 
(VCO) using shift registers and a gate which detects when a certain num- 
ber of pulses are counted. The phase error signal is obtainea using an 
exclusive-or gate whose inputs are the VCO output and a bird-limited 
version of the input. The input to such a digital p.\ase-iocked loop dooc 
not consist of a sampled signal, but rather the hard limited form of the 
FM signal, which provides information about its zero crossings. Other 
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digital techniques [8 3, [9 ] , [10 ]use uniform samples of the FM 
signal but require a voltage controlled oscillator and multiplier tor phase 
detection. The performance of these systems was obtained by computer 
simulation. It is also possible to use a variable sampling frequency 111 I 
to demodulate FM. 

It is interesting to note the digital phase-locked loop progress 
being made abroad. In Japan, hybrid systems have been studied [12 1 
[13 j in which the sampling process occurs after the phase detector, and 
the digital signal converted back into analog form to drive a voltage- 
controlled oscillator. Such phase-locked loops with sampled-data control 
were the forerunners to the all digital phase-locke d loops currently being 
investigated. 

The digital phase-locked loop presented in this thesis operates as 
a (nonlinear) digital filter: the input is the uniform sample sequence of 
the FM signal and the output is the sample sequence of the demodulated 
message. All operations within the loop are digital: gating, storing in 
registers, shifting, adding binary words, and the VOO algorithm. Further- 
more, the system is designed so that it operates on the input sequence in 
real-time, requiring all the digital operations to be performed within the 
sampling period. This feature enables the digital phase-locked loop to 
be constructed and tested using actual modulation, instead of a computer 
simulation. Since the digital phase-locked loop is designed using stan- 
dard logic operations, it can be constructed using LSI, with its advantages 
of low cost and high reliability. 

The structure of the digital phase-locked loop presented in this 
thesis is similar to that of an analog phase-locked loop: the phase error 
between input and VCO is generated by taking the prcduct of these two 
signals, and this error signal is processed by a digital filter. Two dis- 
tinguishing features of this digital phase-locked loop are the real-time 
and synchronous operation. The restriction of real-time operation places 
limitations on the number and type of digital calculations allowable, as 
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well as placing an upper limit on the sampling frequency, as all 
computations must be performed with one sampling oeriod. As a re- 
sult, it is necessary to sample the IF FM signal at a frequency below 
the IF frequency; the permissible sampling frequencies for this sub- 
sampling process are examined in Ch. 2. A second consequence of 
real-time operation is the preclusion of digital ~ | u 1 tipliers, as digital 
multiplication requires excessive computat time Uence ail gains 
appearing in the digital phase-locked loop are powei.. of 1/2, which 
are realized by shifting a binary word. However, we must still construe; 
a phase detector without using a binary multiplication. 

The solution to the phase detector problem is to use a voltage 
controlled oscillator whose wave form is a square wave, having the 
values ±1. Then the binary multiplication of the input by the VCO signal 
reduces to a logic operation: the bits of the input word are exclusive- 
or gated by the VCO output. Since the VCO output is required only at 
the sampling times, an actual oscillator (or counter) is not present in 
the system; instead, an algorithm is developed which determines the 
correct VCO output at the sampling time given all the previous VCO in- 
puts. This algorithm is realized as a logic operation which require? 
negligible computation time. 

Although using a square wave VCO makes real-time operation poss- 
ibile, the square wave introduces undesirable harmonics into the digital 
phase-locked loop. It is impossible to have all of these harmonics fall 
outside the loop bandwidth because ot the frequency aliasing produced oy 
•sampling. In order that these harmonics have a minimum effect on the 
loop operation, a design condition is derived. The harmonic contribution 
to the output is also calculated. 

The digital phase-locked loop internal arithmetic is performed 
with a finite number of bits and therefore quantize iron noise is present. 
An examination of the quantization noise at the error signal leads to u 
second design equation for the loop. 

The original motivation for designing a digital phase-locked loop 
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was to obtain threshold extension over a discriminator. Computer 
analyses of analog phase-locked loops [14 1 predicted up to seven 
dB threshold extension with a third order loop. In order to obtain an 
indication of the digital phase-locked loop performance near -.hrcshoio, 
a deterministic model for an input noise spike is introduced to the loop 
equation, which is solved on a general-purpose digital computer, ana 
the result examined to see whether or not the digital phase-locked loop 
follows the i..put spike. Although this method cannot predict the abso- 
lute location of threshold, it does give a relative perfc mance measure 
between the first, second, and third order digital phase-locked loops. 

The digital phase-locked loops were constructed using DTL and 
TTL logic cards and tested with sinusoidal and constant modulation. 

It is found that the first order loop threshold is identical to thai of a 
discriminator, the second order loop provides threshold extension over 
the first order loop, and, most importantly, the third order ioop does not 
provide any threshold extension beyond the second order loop and it can 
degrade the performance. Tnis important result is found to be caused by 
the third order loop losing lock temporarily in response to an input noise 
spike. 

In addition, to the experimental deter mi . of thresnol i, the effect 

of the number of bits used and the resulting ...'creation on the threshold 
is determined experimentally. It is found that if the sampled input signal 
is coded into a three-bit binary word, threshold is five dB worse than for 
coding into a ten-bit binary word. The result of truncating to ten bits ail 
binary words in the loop is to lose one dB in threshold. 
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CHAPTER 2 

THE ALL DIGITAL PHASE-LOCKED LOOP 

In this chapter, the structure of the digital phase-locked loop 
is developed, with particular attention to the real-time operation 
of the system. The restriction that the unit operate in real-time places 
limitations on the sampling frequency and on the number and type of 
calculations performed by the loop. As a result, a computationally 
simple algor nm is developed to represent a square wave voltage- 
controlled oscillator; but although the algorithm is simple, the 
square wave introduces harmonics into the loop. A condition which 
minimizes the harmonic contribution is developed, and it serves as 
a design equation. A second disign equation involving quantization 
noise is introduced concerning the maximum phase error allowable in 
the digital phase-locked loop. 


2.1 Digital Phase-locked Loop Structure 

The block diagram of the all digital phase-locked loop (DPLL) 
is shown in Fig. 2.1-1. The received noisy FM signal is band-pass^ 
filtered, yielding x(i) , sampled, and converted to a binary word, 

. From this point on, all signals in the DPLL appear as binary 
words. The DPLL error signal, , is the product of the input x^ and 
the digital voltage controlled oscillator (VCO) output w^ . Digital 
filtering of this error signal yields the loop output, y. , ’«hich in turn 
determines the new VCO output, w, , , . The particular dici. al filters 
considered are a proportional path, proportional plus integral paths, 
and proportional plus integral plus double integral paths, yielding 


first, second, and third order DPLL's respectively. The DPLL output 
is converted to a sta.*“case signal by a digital -to analog (D/A) con- 
verter and low pass filtered to produce the analog output, y (t ). 

The digital pha^e -locked loop is d€ signed to operate as a real- 
tima computer, requiring all calculations to be performed In one 
sampling period. This requirement places restrictions on the number 
and type of arithmetic operations possible, and on the sampling frequency. 
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Fig. 2.1-1. The all digital phaso-locked loop. 



2.? Determination of Sampling Frequency 

The digital phase-locked loops were constructed using DTL 
and TTL logic, having speeds of 1 MHz and 10 MHz respectively. 

Allowing for ten logic operations per computation interval (which is 
equal to the sampling period), the sampling frequency must be less 
than 1/(10* msec) = 100 KHz. Therefore it is impossible to sample 
the bandpass signal x(t) (Fig. 2.1-1) at the Nyquist rate or greater. 

It is, however, possible to sample the bandpass signal below the 
Nyquist rate without losing any information. Let x(t) occupy a 
bandwidth B Hz centered about f Q Hz, as shown in Fig. 2.2-la. Then 
it can be shewn [ l4 that an allowable sampling rate is: 

f ^ 2B (2.2-la) 

s 

nf = f ± B/2 (2.2- lb) 

s o 

There are .other possible smapling frequencies. To obtain them, 

first consider x(t) sampled at twice the highest frequency component: 

f = 2(f + B/2). The sampled spectrum appears in Fig. 2.2- lb. If 

the sampling frequency is reduced slightly, spectral overlap occurs. 

But if we continue reducing f g , we reach the situation of Fig. 2.2-lc, 

where f g - 2 (f^ - B/2) and no overlap is present. It is possible 

now to reduce f until the arrangement of Fig. 2.2-ld appears, where 
s 

f = £ + B / 2. Hence, f can fall anywhere within the interval 
so' s 

C(f 0 + B/2). 2(f o -B/2) ] 

Further reduction of f yields overlap, until we reach the situation 

of Fig 2.2- le; here f v f - B/2. We can continue to reduce f until 

so' s 

Fig. 2.2- if results with f = 2 (f +B/2)/3. Therefore, the sampling 

s o 

frequency may fall anywhere in the interval 

[2(f o +B/2)/3. (f o - B/2) 3 

Continuing this process, the next allowable sampling frequency is 

f = 2(f - B/2)/3 snown in Fic- 2.2- lg, and f may be reduced until 
s o # s 

f » 2(f + B/2 )/4 as in Fig 2.2-lh, yielding the interval 
s o 
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Fig. 2 . 2 - 1 . Spectra of the bandpass signal and sampled bandpass signal for various sampling 
frequencies, f . 



[2(f o + B/2)/4, 2(f o - B/2 )/3 ] 

The generalization of this procedure is that the sampling frequency 
must fall in the interval 

[ TTT (f o + B/ 2 >. |( f 0 - B /2>] . (2-2-2) 

with <c an integer, for no spectral overlap to occur. Of course, we 

must also have f ^ 2 B. Notice that the values of f specified oy 
s s 

Eq. (2.2-lb) are included in the endpoints of the intervals of Eq. (2.2-2). 

Summarizing the above, if x(t) is a bandpass signal occupying a 
bandwidth B Hz centered about f Q Hz, then x(t) is specified completely 
by its sample values at instants t = where 

f ^ 2B (2.2-3a) 

f s e [f < f o + B / 2 >- 10 ) u 

[f (f, ♦ B/2). f(t 0 - B/2)] U 

[f(t 0 + B/2). f (£ 0 - B/2)] U 
[srr‘ f o + B / 2 >- r< f 0 - B / 2 >]u 

(21 2-3b) 

The important conclusion is that it is possible to sample x(t) at a 
rate slower than the carrier frequency f Q . 

For ease in implementation, we shall sample according to 

f g * 2Bm (2. 2-4a) 

n.f s - f Q - B/2 (2.2-4b) 

where m and n are positive integers; together, these imply 

f Q /f s * n + 1/4 m (2.2-5) 
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Using m a 1 in Eqs. (2.2-4) leads to an interesting interpre- 
tation. Writing the bandpass signal x(t) as 

x(t) = a(t)cos2nf o t + b(t)sin2TTf o t (2.2-6) 

where a(t) and b(t) are bandlimited to B/2 Hz, the sample values 
are, using Eq. (2.2-5) 

x(kT ) = a(kT )coskn/2 + b(kT )sinkn/2 (2.2-7) 

Note that if k is an even integer, k = 2 p , 

x(2pT ) * (-1 ) P a(2pT ) + 0 (2.2-8) 

s s 

and if k is an odd integer, k = 2p + 1, 

xC(2p+l)T3 = 0 + (-1) P b[(2p + 1)T ] (2.2-9) 

5 S 

Hence the even numbered samples provide information for reconstruc- 

* 

ting a(t), since the Nyquist rate for a(t) is 2 (B/2) = f s / 2, while the 
odd numbered samples provide information for reconstructing b(t). 

Having determined a(t) and b(t), x(t) is determined via Eq. (2.2-6). 

2.3 The Digital Voltage Controlled Oscillator 

A voltage controlled oscillator (VCO) is an oscillator whose frequency 
deviation from its nominal frequency is proportional to the input: 

w(t) = g(2nf Q t + 

where y(t) = VCO input 

w(t) — VCO output 
f * VCO nominal frequency 
g(») = VCO waveform, with g(x + 2 it) = g(x) 

G vco = VCO gain, in (rad/sec)/vok. 

In the digital phase-locked loop we are interested solely in the VCO 
output at the sampling Instants since only these values arc required 
for the computations. The VCO output at the sampling instant t s kT 

3 


Sooj 


y(T) dT 


(2.3-1) 
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is. using Eq. (2.2-5) 


w(k? ) » w. = g(kTT/2m + 


m 

S/CO j 


y(*J) dr) (2.3-2) 


Since we have y(kT g ) a y k available, as opposed to y(t), the 

integration is performed digitally as a summation: 

k-1 

W k * 9 (kTr/2m + 2 y p ) (3.3-3) 

p s - 00 

Mots that the summation extends only up to p = k - 1; this is a 
consequence of the causality of the digital filter. Defining the VOO 
phase by 


5 = ®w», S-„, y > 

we obtain the recursive relation 


®k ~ ®k-l S/OO y k - 1 


and the VOO output is 


k = g(kTT/2m + S fc ) 


( 2.3-4) 


(2.3-5) 


(2.3-6) 


Fig. 2.3-1 illustrates the computations required to generate the 


VOO output. The VCO input y k (identical to the DPLL output) is scaled 
by GyQQ and digitally integrated (or accumulated), forming the VOO 
phase , which is then added to the carrier term, kr/2m. And then 
the VOO waveform g (.) must be evaluated. 

We can avoid the need for calculating the carrier term krr / 2 m by 
recognizing that 

k-1 _ 

W2m + + S Gyoo^ 

p — 0 

k-1 


*0 + f m0 <"/** + G VCO y k ) (2 - 3 ‘ 7) 


That is, the carrier term may be obtained by intc rating the constant 
kTr/ 2 m. Hence, we may modify Fig. 2.3-1 to Fig. 2.3-2, which is com- 
putationally simpler, as now only a constant term is added at each 
calculation. 
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Fig. 2.3-2. The digital VCO 
term, n/2m. 





The next task is to choose the VCO waveform, g ( • ). Two 
factors are important here: the complexity of the g ( • ) computation 
and the multiplication of the VCO output and the DPLL input. 

The function computation must be sufficiently brief as it must be 
performed within the computation interval. This can be accomplished 
using a read only memory (ROM). The function argument is used to 
..address the ROM which is preprogrammed with the function values. 

A typical access time is 22 nsec using TTL. 

The function values determine the complexity of the digital mul- 
tiplication required. For, if the VOO output can assume any value 
(i.e., any of the discrete quantized values), as is the case for g(x) = sinx, 
then the digital multiplier must be able to multiply two arbitrary binary 
numbers, an operation which more computation time than is available. 

One solution to this problem is to utilize a ROM: the binary words to 
be multiplied are used j address the ROM, by concatenating the two 
words, for example, and the ROM stores the product at this address. This 
scheme, however, is made impractical by the memory capacity require z . 
For example, if we are working with 10-bit arithmetic, there are 
2 2 ^ >10 possible products to compute, and hence the ROM 

must store one million 10-bit words. 

A solution to the multiplication problem is to choose the VCO wave- 
form g (• ) to be a square wave, having only values ± 1: 


g(x) = Sq(x) 


+ 1, 0 5x <TT 
- 1, IT <X ^ 2TT 


g(x) = g(x + 2 tt) (2.3-8) 

Hence the binary multiplication is reduced to a simple logic operation: 
if g(x) = + 1 , pass the input with no change; if g (x) = - 1 , form the "one 's 
complement" (or whatever negative arithmetic is used) of the input. This 
operation is accomplished using exclusive-or gates and sufficient compu- 
tation time is available. 
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To complete the digital VCO design, we must specify how g { • ) 
is computed. 

2.4 Digital VCO Algorithm 

Having decided that the VCO outputs are limited to ±1, we must 
now determine which of the two is correct, given the VCO argument, 
k~/ 2m + . That is, we must determine whether the VCO argument 

falls in the interval [0i u ) or Ctt, 2n). Moreover, the intervals are 
considered modulo 2 it to account for the periodicity of g ( • ). But the 
binary numbers group themselves naturally as a result cf the pattern of 
ones and zeros. For example, if we look at a list of the binary numbers 
from 0 to 7, "shown in Table 2.4-1, and focus attention on the second 
digit, we immediately see that this digit partitions the orginal numbers 
into groups of two. Furthermore, we can arbitrarily say that the fiist 
group represents the C 0, tt) interval, the next group represents the 
Ch2tt) interval, the next group represents the [2 tt, S tt ) interval, 
and so on. In fact, if the list of binary numbers is continued above 7 , 
the same identification can be made using the second bit, since this bit 
always exhibits a periodic pattern. Furthermore, if we consider negative 
numbers using offset binary, the same identification can be made, as is 
seen in Table 2.4-2. There is an ambiguity at zero since there are two 
binary representations of zero and hence the VCO output depends on hew 
the value zero is approached. But this difficulty is minor, as the VCO 
argument will rarely be exactly zero. 

Therefore, the VCO output is determined by one bit of the binary 
word representing the VOO argument: if this bit is a 0, *he VCO output 
is +1; if this bit is a 1, the VOO output is - i. Also, to allow for the 
periodicity in Sq(-), the adder used to integrate the VCO input is allow- 
ed to overflow. The particular bit used tc determine the VCO output 
fixes the VOO gain. Above we used the 2's place digit — hence the in- 
terval [ 0, 2 l ) corresponds to the interval L 0, tt ), producing the gain 
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Table 2.4- L The binary numbers and their two's place digit 


Decimal No. 


Offset Binary No. 


2 ; s Place Digit 



- 1 

0110 

- 2 

0 10 1 

- 3 

0 100 

- 4 

00 11 

- 5 

00 10 

- 6 

000 1 

- 7 

0000 


Table 2.4-2. The offset binary numbers and their two's place digit. 
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tt/ 2 1 (rad /sec) /binary number. If the quantization step size is S 
volts* the VCO gain is tt/2 1 S (rad/ sec) /volt. In general, using the 
2 n place digit produces the gain ~f s 2 n (rad/sec)/volt. 1 

Mote" that this gain is implicitly introduced by the \ (. ) calcula- 
tion and so'the-VCO gain is not accurately presented in Fig. 2.3-2. 

Fig. 2.4-1 corrects the situation, where both, an explicit gain, g- ,'as . " 

i 

well as the implicit gain g 2 are shown. The VCG-gain includes both: . 


Gyco = 9i * 92 

Note that we must now add the constant tt/ 2 m with the implicit gain 
removed; hence we add the term (tt/ 2m) /g 2 * Also note that the out- 
put is labeled , the 2 n bit of the integrator binary word. The VCG 

output w, is determined from B via 
k n 


\ s \ 


+ 1, B n = ° 

-1, B n = 1 


If m is a power of 2, the term (7T/2m)/g 2 is simply represented 

4 

example, if the VCO output is determined by the 2 bit (n = 4) then 
TT/g 2 = 10000 
and 

n/ 2g 2 = 01000 


IT 


tt/ 4 g 2 = 00100 
rr/8g 2 = 00010 
TT/lBga* 00001. 


Of course, it is also possible to obtain values which are combina- 
tions of the above, such as 


tr/(8/C)g 2 = TT/2g 2 + Tr/8g 2 = 01010 . 

Finally, in constructing the DPLL the VCO output w^ = ±1 is never 
actually generated. Instead, the B bit is used directly to process the 
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input binary. For if either one’s complement or offset binary arithmetic 

is employed, then negating a binary word is equivalent to complement!,.-, 

it. Hence the B bit can gate each bit of the vat through an exciusive- 
n • 

or gate, as shown in Fig. 2.4-2. If = 0. ga*e cutout is identical 
to the input (a multiplication by + 1 ) while s 3 =1, the gate output 

is the complement of the inpci (a multipli ration by - 1). Hence the mul- 
tiplier output is indeed • x, . Of course, if two's complement 
arithr ' ic were employed, the ex slusive-or gating operation must be 
modified.. 

Summarizing, we see that one bit of the VCO argument serves to de- 
termine the VCO output, and the particular bit used determines the implicit 
VCO gain. The product of DPLL input and VCO output is realized as an 
excldsive-or gating operation, requiring negligible computatic n time . 


■ 2.5 Design Considerations — Maximum Error Signal 

In order for the DPLL to demodulate without distortion, the error 
Signal, e^ in Fir. 2.1-1, must be an accurate measure of the phase error, 
between the Input and VCO. The product of input signal and VCC output 
yields the ternr s’n($ k ~ ^ k ) , and therefore the difference <? k - v k 
must be kept much less than tt/ 2 to have = ‘ < 0 ^ - 0^7 the smaller 
this phase error, the better is the approximation. However, e, r is rapre- 

A. 

sented by a binary word having a fixed number of pits, and decreasing 
its range of values deteriorates the si', sl-tc-quantizafcion noise ratio. 
Hence a trade-off exists between the accuracy or the approxinat_o~ 
sin(9 k - $ k ) ■ a *' c Pj c ~ $ k and the signal-to-quantization noise ratio. 

Let the system A/D converter code analog signals amplitude limited 
to V volts into B-bit binary numbers, producing a quantization step size 


S =- 2 V / 2 


( 2 . 0 -i) 


Let tne phase difference cp^ - range over 2 M of these leve.s: 

I cp k - $ k [ S MS y f! . s >— 2 1 
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Then assuming v, - is uniformly distributed over its range, tht 

K ii 

quantization. noise is 


N = S" / 12 

<4 

and the signal power is 

- ,2 - i ^ , 2 ‘ :.i<WDqy*+i) 3 r 

E ^k " ~ k 5 zx-r ; " v (n ~' " 3 


v ~.0-o) 


(2.5-4) 


The signal -to-euantir ation noise ratio is then 


3 / N_ = 4 M { M * i ) 


(2.5-5) 


As a re suit of the aoproxinaticn sir. (o _ - © } °* c - c- the 

v k k' k k 

harmonic noise introduced is 
N d ■ 


M 2 

TmTi ~ (nS-sinnS) 
n=-M 


and the signal -to-harmonic noise is 


(2.5-S) 


_S_ = MfM+IH2M+l)S 

N , M 

3 £ (nS-sinnS) 

n= - M 


( ft Z ^ ' 


Eqs. (2.5-5) ar.d (2.5-7) are plotted in Figs. 2.5-1 ana 2.5-2 as a 
function of the maximum phase difference, M S, for the case of an A/D 
converter which accepts =5 volts maximum end uses 1C and 12 bits 
respectively. As expected, the signal -to-quantization noise ratio improves 
when the signal range increases, while the signal -to-harmonic noise re.io 
deteriorates for increasing signal range. The die -.tat phase-locked loops 
are designed so that these two ratios are equs' . From Figs. 2.5-1 ar.d 2.5-2 
the maximum allowable phase errors are: 
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Fig. 2.5-1. signal -to- quantization noi.'tc nnc! Htynal -to -harmonic noise oa a function of the maximum 
phasf; cllf foremen for n ton-bit A/D convertor, 



Nk mber of Bits 

Maximum Phase Error 

10 

0.35 

12 

0.23 


Table 2.5-1 


The maximum allowable phase error places a restriction on the loop 
gains as will be seen in the following chapter. 

2.6 Design ^cnstdearattons— VOQ Harmonies 

Although tire square wave function of Ec *--3-8) adopted for the 
VGO leads to a simple compuretianal algorithm, it has tire disadvantage 
of introducing^ harmonics into the DPLL. For the fmctlon Sq(x) has the 
Fourier series 

Sq(x) * ^£sinx + ~sin3x + ~sin5x +'••• + ~j^sin(2n+l)xi-- j (2 

and the product -2ooe C2uf 0 t + P(t)]Sq[2nf o t+#(t)l generates 

Imwaiicr at f * 0, 2f , 4f . 6f . .... 2pf Now in an analog 

o o o o 

phase-locked loop* these harmonics would be sufficiently attenuated 
by tiie loop filter and tin VOC integratioa. However, in the digital phase 
Indeed loop* the samgjHag process shifts these frequencies, and eventually 
a harmonic will be folded (aliased) into the loop bandwidth. Since 
f Q « (n + l/4m)f s , the frequency f Q gets aliased to f s /4a, and therefore 
" the iradtjpUer output contains harmonics at 

0* f /2m, 2f /2 m* .... pf /2m* ... Hz. 

8 S 8* 

th 

But the frequency f is equivalent to zero frequency and so the 2m 

3 

harmonic is aliased to d*c.* falling within the loop bandwidth. The pro- 
d uc ts 

[-2 cos (kn/2m + cp^ ) J ^£~ 7 yj-sin [k (4m-l)" /2m + (4m-l)$ k ] v 2 

and 
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[-2 cos (Jen /2m ♦ <P fe ) ] ~ ~rj sin [k(4m+l)^2m + (4m+l)$ k ] (2.6-2b) 

Ak 

generate this 2ia“ harmonic which is 

o {l^l C (««-»)» k * »,J ♦ sin C (4aw-I)» k - * fc } } (2.6-3) 

If tire loop is following the input modulation, then $ . w ¥, and the 

th k K 

amplitude of this 2m harmonic is approximately 

If -JL * -JL } . S|m 

She afove analysis Indicates foal m should be made large in order 
that the aliased harmonic base small amplitude and not interfere with 
the loops tre&dhg ability. How e v er, the multiplier output first harmonic 
is at f » f / 2 m. and increasing' m causes ibis harmonic to fall within 
the DPLL bandwidth. As tire amplitude of the first harmonic is approxi- 
mately 4 (1* l/3)/tr ■ 16/2ir. it is tire first harmonic that poses the 
- V' -- - -1 * th 

greater threat of loop malfoaction tiren tire 2m harmonic. If the input 
modulation produces a frequency deviation Af Hz. foen tills first harmon- 
fre< * lancy (assumtaa \ « * k > 2if 

Hz. Requiring that this first hajnmonic foil outside the loop bandwidth. 

.we have • 

(f e /2m) - 2Af (2.6-4) 

Eq. (2.6-4) and Table 2.5-1 form tire basis of design of the first 
order phase-locked loop. Chapter 3. Consideration of the DPLL* s re- 
sponse to input spfkes tireh allows design of the second and third order 
DPLL's, in Chapters 4 and 5. 


2.7 






'While the use of a square wave VCO simplifies the DPLL implementa- 
tion. it introduces a severe nonlinearity into the DPLL difference equation, 
making it impossible to solve exactly exoept In a limited number of 
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special cases (the first order DPLL equation is solved in Ch. 3 when the 
Input carrier is unmodulated).. We therefore introduce a linearized model, 
developed on the basis of tire design conditions of Sections 2.5 and 2.6. 

The product of input carrier and square wave VOO generates harmonics 
at O', f /2m, 2f /2m, ...» pf /2m, ... Hz, with p an integer. Of 
these harmonics, only the d.c. term contains information about the phase 
error between input and VOO, and ideally, fids shouH be the only term 
present. Eq. (2.6-4) stipulated that the additional harmonics fall outside 
the DPLL bandwidth. As previously ^observed, the sampling process folds 

Aft. 41, 

the 2m“* harmonic to dw04 also, this 2ni harmonic has an amplitude 
of approximately &m/(16m* - 1) (relative to tee d.c. term). Hence, to 
obtain a linearized model we may neglect all the harmonics except the 
first, so teat 

\ - £«*«(** -* k > . Pw-U 

Furthermore, tee phase difference (<*^ - m fc ) is restricted according to 
Table 2.6-1, and we there fo re can further approximate: 

•k “ nt*k**k> < 2 ‘ 7 ' 2 > 

Using Eq. (2.7-2), the linearized model of tee DPLL appears in Fig. 2.7-1. 
Note that, tee input to tee linearized model is tee input phase, q>^. 

The suppressing of the harmonics is supported by the square wave 
nature of tee VOO: if tee VOO argument lies anywhere in tee C 0, tt) 
interval, the VOO output is a * 1. Therefore, although the true VOO phase 
is perturbed from teat in the linearized model, tee VOO output sequence may 
still be identical. 

2.8 Hardware Considerations 

The digital phase-locked loops designed were constructed using 
commercially available DTL and TTL logic cards. As previously men- 
tioned, all computations me performed within one sampling period, which 
is chosen to be 20 psec to accomodate tee logic speed. All computations 
are performed in parallel form. 
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The A/D converter employed generates words in offset binary 
form: the most significant bit is the sign bit which is a 1 for positive 
voltages; negative voltages are represented by the complement of the 
corresponding positive voltage word. There are two representations of 
zero volts: +0 * 100...0; -0 *= 011...1. Table 2. 8-1 illustrates the 
offset binary words obtained when the input voltage is between -7 and 
+7 volts: for simplicity only 4 bits are used. The A/D converter used 
in the DP LL converts voltages between -5 and +5 volts into 10-bit 
offset binary words. 


2.8.1 Multiplication by 1/2** 

To avoid a binary multiplication (which is excessive in computation 
time) all gains appearing in the DPLL are chosen to be powers of 1/2. 

As a result, the multiplication operation becomes a shifting operation, 
shifting the binary word once for each factor of 1/2. However, care must 
be taken in performing this shifting operation: The sign bit must NOT be 
shifted, and the bits which become vacant mast be tilled with the comple- 
ment of the sign bit. For.example, 


4 (+*) * 4 (U00> « 1001.00 

and 

£ (-4) * 4 ( 0 °!!) * 0110.11 


(if the answer is allotted only four bits, then the two bits to the right of 
the binary point are dropped. We shall see later that in the DPLL' s 
designed here these bits are net dropped from the answer.) 


2.8.2 Addition Algorithm 

In order to correctly add two offset binary words, an algorithm is 
required. Let the two N-bit words to be added be denoted by 
A = a i 2 . . . a M and B ■ b t . b^, with a x and bj the 

sign bit. The algorithm consists of the following steps: 
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Voltage Levels (volts) 


Offset Binary Word 


+ 7 
6 
5 
4 
3 
2 
1 

0 

- 1 
- 2 

- 3 

- 4 

- 5 


6 


- 7 


{ 


1111 

1110 

1101 

1100 

1011 

1010 

2001 

1000 

0111 

0110 

0101 

0100 

0011 

0010 

0001 

0000 


Table 2.8-1. Illustrating offset Unary coding. 
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1* ?orm tlio ordinary binary sum A+B = c o c l c 2*** c N* where 
c . = the sign bit cany. 

i* 

2. Use the carry into the sign bit as an ervd-around carry, resulting 

in the new sum d 0 d i d 2*** d N* 

3. The answer in offset binary is obtained by deleting the digit d ^ : 

A+B ■ d 0 d 2 d 3'--V 

Several examples will elucidate the algorithm: 

Example 2 .8-1 +3 * 1011 

+2 * 1010 

Step 1. Q ^ carry into sign bit 

1011 

1010 

10101 

Step 2. 10101 

0 

10101 

Step 3. Answer = 1101 = +5 

Example 2.8-2 +3 * 1011 

-2 - 0101 

Step 1. j^-carry into sign bit 

1011 

0101 

10000 

Step 2. 100 00 

1 

10001 

Step 3. Answer “ 1001 * +1 
2.8.3 Arithmetic Saturation 

If two numbers whose sum exceeds the voltage range are added using 
the above algorithm, an erroneous result occurs, referred to as overflow. 
For example, + 4 + 6 * 1100 + 1110 * 1011 (using the algorithm) 
but 1011 - +3. 
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In linear digital filters it is possible for overflow to generate limit 
cycles, which can be eliminated by having the addition process saturate C 16 J 
The DPLL is a nonlinear digital filter, and to prevent the possibility of 
limit cycles, the arithmetic is modified to include positive and negative 
saturation when the r um exceeds the allowable limit. 

First, we must detect when overflow occurs, and whether it is positive 
or negative. Note that overflow can never occur when adding numbers of 
opposite sign. When two positive numbers are added, the sign bit carry is 
1; if in addition, overflow occurs, the carry into the sign bit is also a 1. 
Hence positive overflow occurs when both of these carries are 1. Adding 
two negative numbers always produces a 0 sign bit cany, and overflow 
produces a 0 cany into the sign bit; hence negative overflow occurs when 
- both these carries are 0. These two conditions are easily detected using 
nand gates, as shown in Fig. 2.8-1. 

Having determined an overflow situation, we may generate the most 
positive or most negative number, whichever is required. The adder plus 
saturation processing is shown in Fig. 2.8-2, where the adder is followed 
by a storage register and a bank of nand gates. The inverted output of the 
storage register is used; hence R = the complement of whatever is stored. 

The negative overflow signal is inverted and used to reset the register; the 
positive overflow signal gates every bit through a nand gate. 

When no overflow occurs, N * P * 1, and the register is not reset. 

The register stores C and its output is C; each bit is complemented by 

ear 

the nand gates producing Y ** C « C. 

When positive overflow occurs, P 3 0 (N 3 1) and the nand gate 
output is Y * 1111, the maximum positive voltage, regardles s of the 
register output. 

When negative overflow occurs, N » 0, P = 1 and the register is 
reset to zero, producing R » 1111. The nand gates invert this producing 
Y * 0 0 00, the most negative voltage. 
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Fig. 2.8-1. The required logic for detecting positive and negative overflow. 
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2.8.4 VOO Overflow 


The VCO algorithm was developed in Sec. 2.4 where it was concluded 
that the VCO integrator must be allowed to overflow in order to realize the 
periodicity of the square wave function, Sq ( * ). This means that if, for 
example, we are working with ±7 volts maximum*., the addition must be per- 
formed on a modulo 8 basis. 

Let us examine the addition algorithm of Sec 2.8.2 by adding +5 and 
+4. The result should be 

*5+4 = +1 mod 8 (2.8-1) 

The addition algorithm, without saturation of course, yields 

1 

+5 * 1101 
-*-4 = 1100- 
11001 
* m 

11010 

+ 5+4 » 1 1010 = +2 

Note that the answer is too large. In fact, the addition algorithm will 
always generate a larger result because the overflow (the cany into the 
sign bit) Is the end around carry. Hence to achieve the correct overflow 
sum, we must prevent the end around carry operation. 

Next, consider the sum - 5 - 4, whose sum is 

- 5 - 4 = - 1 mod 8 

The addition algorithm yields 

0 

-5 * 0010 
-4 » 0011 
00101 

0— J 

00101 

-5-4 * 0101 = -2 

Here, the difficulty is that the end around carry is a zero, when it should be 
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a one. This will always be the case for negative overflow because the 
operation -A - 3 is really the sum (111! - A) + (1111 - B) 

■— ililO-A-B, ana for negative overflow the end around carry is zero, 
so that our result is different from the correct result, lllli-A-B. 

Therefore, to obtain a correct result for the VOO integration (addition), 
it is necessary to make two modifications in the existing addition process. 
First, the saturation hardware must be eliminated. Second, the end around 
cany must be complemented whenever overflow is detected. The end around 
carry modification may be simplified by noting that under normal loop opera- 
tion the VCO integrator will never overflow negatively, since the positive 
number 'n/ 2 m is added at each computation. Hence, the end around cany 
may be set to zero for all VOO computations. 
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CHAPTER 3 


THE F.. ST ORDER DIGITAL PHASE-LOCKED LOO? 


3,1 Tha First Order DPLL Equation 

The first order digital phase-locked loop is shown in Fig. 3.1-1. 
Referring to this figure, let the sampled FM signal plus noise be de- 
noted by x^: 

^ * -2cos(2Trf Q kT s + c^) + n^ (3.1-1) 

where 

f Q = carrier frequency 

T = sampling period 
s 

* input phase 
n^ * IF noise 

The reason for choosing the signal amplitude of 2 volts is seen below. 
Incorporating Eq. (2.2-5), we have 

» -2cos(kn/2m + + n^ (3.1-2) 

The error signal e^ is the product of the input and VCO output, w^: 

«k " v w k (3a - 3> 

The DPLL output, y^, is simply proportional to the error for a first 
order loop: 


y k * 9 ' e k 

where g * forward loop gain. The VOO phase is then 


(3.1-4) 


\-l + ‘Voo^-I < 3 - l - 5 > 

where * VOO gain, which includes both any explicit scaling 

of the loop output, as well as the Implicit gain introduced by the Sq (• ) 
algorithm, as explained in Ch. 2. Finally, the VCO output is 


w - Sq(kn/2m + ) 


(3.1-6) 
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Fig. 3.1-1. Block 





Combining Eqs. (3.1-2) through (3.1-6), we obtain the equation for 
the VCO phase : 

l * © k - 2 G cos (kTT/2m + ^) Sq (^210 + 6^) 

+ Gn^Sq (kTt/2m + ) (3.1-7) 

where G s G^^g * loop gain. To bring this equation into a clearer 
form, expand the Sq ( • j function using a Fourier series 


Sq(x) * - 2 —7 sin ( 2p+ 1 ) x 

p»0 


(3.1-8) 


and expand the product of the fundamental component and the input 
FM signal. Then 

-2Gcos(kn/2m+ <P k )Sq(kr/2m + 6^) « sin (o^ - $ fc ) 

+ ^Gsin(kr^2m + 9^ + ^) 

CD 

-2Gcos(kTT/2xn + <P.)^ £ sin(2p+l )(k^/2m+^) 

p«l ^ 

- losing -» k >* £ G W*k> 


(3.1-9) 


The first term is a measure of the discrepancy between the input and 
VCO phases while h^ is the contribution from the harmonics generated 
by the square wave VCO. The DPLL equation is thus 
4 


$ k+1 -$ k + £ Gsin($ k - © k ) * Gn Jc Sq(kTt/2m + $ k ) 

* ? G WV 

The left side of this equation is the digitized version of a first 
order analog loop, where the. derivative is replaced by a difference: 

df * 1 + G «uao a * taC4 «-' ,, < t > 5 

~ < Vl - *k >/*. - G «Ml 08 ,in<S k - V 


(3.1-10) 


(3.1-11) 
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Ncia zhc.t the corresponding analog loop gain is 


G 


analog 


n f G 
TT s 


<3.1-12) 


The factor 4 /tt represents the Fourier coefficient of the square wave 

fundamental while the factor f enters because computations in the 

s 

digital loop arc performed every sampling period, T . 

s 

The first order DPLL equation, Eq. (3.1-10), is a first order, non- 
linear difference equation, whose solution is obtained under special 
input conditions in Secs. 3.3 and 3.4. 


3.2 Stability Condition or the Loop Gain 

A first order analog loop possesses stable VGO states (phases) 
for any positive loop gain when the input carrier is unmodulated 
(o^s 0) and no noise is present 0). For any initial VCO phase, 

the VOO phase eventually settles to zero (modulo 2it); the larger the loop 
gain, the smaller the transient time required. The digital loop has differ- 
ent restrictions. 

First, consider the case shere the added harmonics h^ are not in- 
cluded, yielding the equation 

Vi = - * G8ln ®k (3 - 2 - l > 

In order that — - 0 for any initial value P q , we must have 
( 4 /tt) G < 2, as demonstrated in App. 1. larger loop gains require smaller 
transient times (where transient time is defined, for example, as the time 
required for the VOO phase to be within 0.1 radian). 

Now if the harmonics are included, the gain restriction changes, as it 
is no longer possible to have 0 In the steady state. For example, 

assume that P^ m 0; then 

p fc+1 - ^ * G e^ * -2G( cos kn/2m ) ( Sq kn/2m ) (2 .2-2) 

which is shown in Fig. 3.2-1 for the case m « 4, and when this integrated 
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Fig. 3.2-1. The assumed DPLL error signal e^ when both input and VOO 
phase arc. identically zero. 



Fig. 3.2-2. The steady state error signal generated when the input 

» 

carrier is unmodulated. 
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to form ^ , the result is not identically zero, as originally assumed. 

Therefore, more care must be taken in calculating the VCO phase. 

Assume that at k = 0, the VCO phase is zero: ® =0. Then 

• o 

e = - 2 cos 0*n/2m Sq(0-:/2m + 0) = -2-1 = -2 
0 


and 

®, = ® + Ge * 0 - 2G 

1 o o 

Now the VCO output should be + 1: 

w 1 = Sq(l*ir/2m+ = Sq (n/2m - 2G) = + 1 
Therefore 

n/2m - 2G ^0 


and 

G — ^/4n (3.2-3) 

We shall see that this (in addition to the obvious requirement G > 0) is 
precisely the gain restriction for stability. Then 

ej * (-2 cos -/ 2m) • (+ 1 ) 


and 


q> 2 = 0 - 2G - 2 Geos tr/ 2m 


implying 

w 2 = Sq(2n/2m - 2G - 2Gcosrr/2m ) = +1 


since 

2n/2m > 2n/2m - 2G - 2Gcosn/2m = (TT/2m - 2G)+(TT/2m- 2GcosTT/2m) >0 
In fact, it follows that for 0 ^ k ^ 2m, the VCO output is + 1, since 


k- i 

£ = <p + G 2 -2cospTr/2.s 

k ° p=0 

and 


0 ^ktT/ 2m + < tt for 0 £ k^ 2m 
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At k = 2r. t 1 , the V CO phase is 

2 m 

v 7 . r= s? + G 23 — 2 G cos p” /2m = 0+0=0 

~ v "‘ ° p-0 

and hence the VGO output is - 1: 

W 2m+1 = Sq((2m+l)n/2n + 0) = -1 

Similarly, the VCO output remains at - 1 for 2m+l ^ k ^ 4 m- 1 . Finally, 
at k = 4m, the VCO phase is 

2m 4m-l 

+ G £ -2cospn/2m + G £ (-l)(-2cosp~/2rr. } = 0 

p = 0 p=2m+l 

ar.d the sequence repeats. The resisting steady state output is shown in 
Fig .3. 2-2 for the case m = 4; it is an eight-point sequence followed by a 
seven-point sequence, etc. Note that this signal has a zero average 
value, while the sequence of Fig 3.2-1 does not. 

The gain restriction, Eq. (3.1-4) has another interpretation: t pre- 

vents the VCO output from executing two consecutive sign changes; i.a., 
it prevents VCO jitter. To illustrate this, note that the VCO argument al- 
ways increases: 

A(kTr/2m + cp fc ) = TT/2m + A = n/2m tGe^ > ~/2m - 2G ^ 0 

Therefore, if the VCO argument crosses from a - 1 to a + 1 interval, it 
must remain in the + 1 interval at least for the next sample. In fact, the 
most the argument can increase is Tt/ 2 m + rr/2m = tt/bi, aud hence the 
VCO output remains + 1 for at least m samples. 

Fig. 3.2-3 shows photographs of the D/A outr ■ : ;he first order 

DPLL for different loop gains. The DPLL has un .c'.vurd loop gain 
(g = 1 ) and an explicit VCO gair of 1 / 16; th ...... cit VCO gain was 

varied, producing different loop gains. The A/D converter codes a = 5 
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(C) 6 s TC / 5 

Fig. 3.2-3. D/ A outputs of a first order 0PLL having loop gain G and 
an unmodulated input carrier. 
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volt input signal into 10-bit binary words; the sampling frequency is 
50 Klin and m 38 4. 

I'ig. 3.2-3a shows the error signal when the VCO output is chosen 
using the third most significant bit of the VCO argument. Since this 
bit carries the value 5/4 = 1.25 volts, the implicit VCO gain is tr/i.25 
and the loop gain is 

G = (1/16) (tt/ 1.25) = n/20 

Note that this gain is less than the critical gain, Tt/4m = tt/16, and so 
the DPLL should generate the sequence of Fig 3.2-2. Indeed, examining 
Fig 3.2-3a we see exactly the anticipated sequence; nine levels, follow- 
ed by seven levels, followed by nine levels, etc. The spikes seen on 
the left side of some of the levels is a flaw in the D/A converter, not in 
the DPLL. 

% 

Next, in Fig 3.2-3b, the VCO output is chosen using the fourth most 
significant bit, yielding an implicit gain of tr/G. 625 and a loop gain 

G = tt/ 10 

which is greater than the critical gain and therefore should result in an 
unstable DPLL. Fig. 3.2~3b bears this out: after two negative levels, 
the VCO output changes sign, producing a positive level, and changes 
sign once again, generating a negative level. This is precisely the con- 
secutive VOO sign changes referred to above. Such a DPLL does not gen- 
erate the desired steady state sequence and is classified as unstable. 

Finally, the loop gain is again doubled by using the fifth most signif- 
icant bit for the VCO output; G - tt/5. Fig 3.2-3c shows that now there 
are two occasions where the VCO output changes twice consecutively; the 
third and eighth levels are inverted and the VCO cc -sets itself on the 
fourth and ninth levels. 

One final point concerning the steady sta. Crignally it was assumed 
that the initial VOO phase, , was zero, which led to the steady state 
error sequence. But does not necessarily have to be zero to generate 
this sequence. In order that the initial VOO output be + 1, we require 
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0 £ V £ it 
o 

To assure that the VCO generates 2 m successive +1 outputs, we examine 
its argument at k * 2m; if this argument is less than then all VCO out- 
puts for 0 — k £ 2m are indeed +1, since the VCO argument always increa- 
ses. Hence we require 


But 


2mn/2m * ^ 


2m-l 

$ s $ +G 2 - 2 cos pfr / 2rr. * $ - 2G 

2m o „ 1 o 

p=0 


and therefore 


Therefore, if 


V < 2G 
o 

0 S $ £ 2G 

o 


(3.2-4) 


we are in the steady state. 


3.3 Transient Response to an Unmodulated Carrie : 

An analog phase-locked loop having positive loop gain will drive any 
initial VCO phase to zero (modulo 2 tt) in response to an unmodulated 
carrier input. Initial VOO phases of ± tt (modulo 2 rr ) place tne loop at an 
unstable equilibrium point and theoretically the loop should remain there 
indefinitely. However, in practice, it is impossible to remain preciesly at 
this point and once perturbed from it * however slightly, the loop migrates to 
a stable equilibrium point. Theoretically it requires infinite time for the 
VCO phase to reach zero; in practice, of course, this transient time is finite. 

The VCO phase of tire digital phase-locked loop desc-ibed by Eq, (3.2-1) 
(i.e., where the error signal harmonics are suppressed) is also driven to 
zero for any initial VCO phase, provided the gain restriction 0 <40/^ <2 
is met. Theoretically, infinite time is required for the VCO phase to roach 
zero, except for certain isolated initial VCO phases, for which a finite time 
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is required. This is pursued further in App. l. 

V/hen the square wave VOO harmonics are included yielding Ec. (3.1-7) 

the D?LL reaches the steady state output sequence for any initial VOO 

phase (including n modulo 2 tt}. Furthermore, the transient time is always 

finite, depending on the loop gain and the initial VCO phase. 

V/e have previously seen that if 0^6 < 2G t..ere is no transient — 

o 

the steady state sequence begins immediately. If > 2G the VCO out- 
put changes sign before k = 2m + 1, widen reduces the VOO phase in a*i 
effort to reach the steady state. More precisely, if 

2G <P <2G + 2GcosTr/2m + Ti/2m (3.3-1) 

o 

then the VCO output w is 

C +1, 0 — k — 2m-5* 1 

W k " ^ - 1, k = 2m (3.3-2) 

This result is proved in A?p. 2. The VCO phase at k = 2m is 

2m-l 

= 9 q + 6 2 w -.( _ 2 cos prr/2m) 

P- 0 

2m- 1 

= q> +G S (+ 1} ( - 2 cos pTT/2m ) 

0 p=0 

— cp — 2G (3.3—3) 

o 

That is, the VOO phase is reduced by 2G. Therefore, if in addition to 

Eq. (3.3-U we have $ <4G, then <¥> < 2G, and steady state is reached 

o 2m * ^ 

after 2m iterations. If 4G — < 6G, then 2G ^ anc * 

0 < p < 2G, and steady state is reached after 4m iterations. Continu- 
4m 

ing in this manner, the region in the - G plane specified by .3.3-1) 
is partitioned into smaller regions? each region has a finite transient time 
associated with It. Fig. 3.3-1 illustrates these regions for m = 4. 

If 

2G + 2GcosTr/2m + Tr/2m ^ < 2G+ 2Gcos -/2m + 2Gcos 2Tr/2m+ 2Tr/2m (3.3-4) 
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f C < k ^ 2m - 2 

W k = j - 1, 2m - 1 < k< 2n (3.3-5) 

2); i.e., the VCO output changes sign at k = 2m- 1. The VCO phase 
at k = 2a is 


2ra-2 

$ * $ + G 5 (+1) (-2 cos p-/2m) + G (-1) (-2 cos (2a-l) n/2ia) 

° p=0 

** - 2G - 4 Geos ^/2m (3.3-6) 

Note jan Eqs. (3.3-4) end (3.3-6) imply that 
- 2Gcosrr/2m + n/2a S p 

and since 2G < w/2m ( the stability condition ), 

0 <5 2m 

so that we must fall into a region already considered. In fact, we can 
partition the region of 2q. (3.3-4) by shifting the previously established 
-regions upward by the amount 2G +• 4GcosT?/2m, and then shifting the 
newly created regions until the entire region, of Eq. (3.3-4) is partitioned. 

The transient ome associated with a shifted region is simply 2m 
greater than the transient time associated with the original region. The 
original and shifted regions are illustrated in Fig. 3.3-2. 

Generalizing, if 


r-1 r 

2GcospTt/2m + (t-1)tt/ 2© S $ <£ 2G cos p^m + r^/2m (3.3-7) 

p=0 ° p--C 


then 

( +1, 0 <k ^ 2m - - 
W k * ] -1, 2m-r <k - 2m 

where OSr ^ 2m (App. 2). The VCO phase at k-- 


2m is 


(3.3-8) 
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(-l)(-2cos ptT/2m) 


2m 


2n-r 

6 +G £ (+l)(-2cos pTt/2m) 

° p=0 


273—1 

+ G £ 

p=2m-rs-l 


r-1 

= $ - 2G - 4G £ cos pn/2m 
° P=1 


Again, 


0 < 



< 6 

o 


(3.3-9) 


(3.3-10) 


(App. 2) and therefore the transient requires a finite number of iterations. 

The duration of the transient for a given and G is determined by 
building upon the regions of Fig. 3.3-2 until the entire rectangular re- 
gion 6 < G < n/4m, 0 ^ <P ^ n is partitioned. For m * 4, the result 

o 

of such a construction is illustrated in Fig. 3.3-3. 

Next, we consider -tt S $ <0. This will complete the analysis 

o 

of the transient response since the DPLL responds identically to as it 

does to $ + 2^n, with n an integer. When $ is negative, the VGO out- 

o o 

put is initially - 1, which in turn produces a positive output which UICTcEScS 

the Initially negative VCO phase. However, unlike the case for positive 

$ , it is possible to fall into the steady s J ate before 2 m iterations. Before 
o 

beginning the transient analysis, we list the VCO phase in the steady state 

secuence- This list appears in Table 3.3-1. Recognizing that the steady 

state sequence may begin at k = 0 or k = 2m, the steady state VC C Phases 

intervals become the union of the intervals for Q and $ , 0 , listed in 

n n+2m 

Table 3.3-2. Hence, if at any time during the transient the VCO phase falls 
into the interval given in Table 3.3-2, then the steady state sequence is 
generated. 

U 

-2G - ir/2m S cp Q <0 (3.3-11) 


then 


* 




< 26 


0 £ 9 

o 

-26 * < 0 


n-1 

- 26 S cos pTT/2m S ^ 

p-0 n 


n-1 

< 26 - 26 S cos pn/2m 

p=0 


■*2G 5 9 


2m 


<? 0 


0 * x *e 


-2Gcosn/2m ^ _$ 


2m+2 


< 2G - 2 Geos Tr/2m 


n-1 " ' ; n-1 

r2GS cosp*r/2m ^ - < 2G - 2G £ cos prc/2m 

p=l , n p=i 



0 £ $ < 26 
. 4m 


Table 3.3-1 The VOO phase in the steady state D?!L output sequence. 


SO 



-2G £ 9 < 2G 

o 

-2G ^ $ < 2G 


n-1 

2G 2 cos pn/2n: £ ^ 


n-1 

< 2G - 2G £ cospn/2xa 
p“l 


_ -2G ^ < 2G 

221 


Table 3.3-2. The intervals that result from taking 

interval and the $ , _ interval of Table 3.3-'.. 
n+ 2m 


cr.e union of the 9 

n 


SI 



(3.3-12) 



-1. k = 0 

- 

+ r. l S k ^ 2m 


(proved in App.2). Thai is. the first VOO output is - 1, causing the VOO 
phase to increase, and the VOO to generate a + 1 output at the next 
sample time. Also, 


2m- 1 

$ >9 + (-l)(-2cos0 ir/2m) +G £ (+l)(-2oosp*?/2m) 

2m O psrl 

s$ + 26 (3.3-13) 

o . - - , 

and therefore 

* -26 « -w/2m S < 26 


ann so we have certainly reached steady state after 2m iterations. How- 
ever, we actually reach steady state sooner; for if -26 — < ®. 

T able 3.3-2 tells us that we are la the steady state. Furthermore, since 

m $o +■ 6 (-l)(^2cos0tr/2m) « 26 

we re ac h the steady state after one iteration when 

-46^P a < -26 

Continuing, if -66 £ < -46. the after 2m iterations. 

-46S^ 2m <-26 

and hence the total transient requires 2m+ 1 iterations. If - 86 -$ Q < - 66, 
the transient requires 2m + 2m + 1 ** 4m + 1 iterations. Ext en di n g these 
results, the region of Bq. (3.3-11) is partitioned as shown in Fig. 3.3-4, 
where each region is labeled with the associated transient time and m « 4. 

If ' 

- 2tt/2m - 2G - 26dosn/2m ^ < -2G ~tT/2m (3.3-14) 
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osksi 


than 

k +1, 2 ^ k £ 2 a (3.3-15) 

(App. 2}. Hare the first two VCO outputs are - 1 before the VCO out- 
put changes sign. The VCO phase at k « 2a is 

1 2a- 1 

$ « ^ +GS (-l)(-2cospTr/2a) + G S (+l)(-2cospn/2») 

p=0 p=2 

« + 2G + 4GcosTr/2m (3.3-16) 

and combining this result with Eq. (3.3-14) we obtain 

$ £ < 2G (3.3-17) 

so that steady state is indeed reached. Here it is possible to fall into 
file steady state sequence at k = 2: from Table 3.3-2, this happens when 

-2G - 2GoosTT/2m £ ^ < 2G - 2Gcos^/2a (3.3-18) 

Since « $ + 25 + 2Gcos^/2a, Eq (3.3-18) becomes 

-4G - 4 Geos rr/2m ^ 6^ < -4Gcos^/2m (S.3-19) 

The intersection of the regions of Eqs. (3.3-14) and (3.3-19) yields a 
region which requires two iterations to reach steady state. The remaining 
portion of the region of Eq. (3.3-14) is partitioned by shifting the regions 
of Fig. 3.3-4 down by 2G + 4Gcosn/2m, as specified by Eq. (3.3-16); the 
result is illustrated in Fig* 3.3-5, again for the case m - 4. 

The generalization of the above is that for 


- F(r+1) ^/2m + £ 2G cos ptT /2m 1 

P*0 


p=0 


(3.3-20) 


and <P ^ -tt, the VCO output is 
o 
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-1. 0 s k < r 

+ 1, r + 1 ^ k ^ 2m (3.3-21) 

where 0 £ r ^ 2m - 1 (App. 2). The VOO phase at k * 2m is 


r 

^2m S ^o + G S ("W- 200 ^ 11 / 2 ®) 

P=0 . 

2m- 1 

+ G 2 (+l)(-2cospTr/2m) 

p=r+l 


= $ + 2G + 4G 2 cosprr/2m 

° P=1 


<p < & < 2G 

o 2m 

(App. 2) so that the steady state is approached without an overshoot into 
the $> 2G region. Steady state is reached at k = r+ 1 if, by Table 3.3-2, 

r r 

-2G 2 cospTT/2m ^ $ <26 -20 2 cospTi/2m 

p=0 r 1 p=l 

Since 

r 

$ r+1 B + G 2 ( -1 )( - 2 cos pTT/2m ) 

p=0 

Eg. (3.3-24) can be expressed in terms of : 

r r 

-4G 2 cos pTT/2ra ^ $ <4G-4G2 cospn/2m (3.3-26) 

p=0 ^ p=0 

Note that it is impossible to fall into the steady state at k = r+ 1 if 
r > m, since the regions of Eqs. (3.3-26) and (3.3-2G) do not overlap: 


(313-24) 


(3.3-25) 


(3.3-22) 


(3.3-23) 
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-4G S cos?~/2m > - ; r~/2m + £ 2Gco$prr/2m ! for : > m 

p=0 c p=0 

Vv'a may now proceed to construct regions in the - G plane and 
label the required transient tine. We start by confining attention to a 
region specified by Eq. (3.3-20) for a particular value of r; let us call 
this region R. The. first step is to see if it is possible to fall into the 
steady state sequence after r iterations; i.s., we find the intersectio. 
of the regions of Eq. (3.3-20) and R. (As noted previously, this intersec- 
tion is null if r > m.) This region requires r transient iterations. For 
any other point in R, the VCO phase satisfies Eq. (3.3-22); i.e., after 2m 
samples, it increases by the amount 

h = 2G + 4G 2 cosp~/2m. 

P=1 

Hence, we translate the previously constructed regions upward by this 
amount, consider only their intersection with R, and observe that the 
transient time for the translated region is 2m greater than that for the 
original region. After this construction, there will still be a portion of 
R which is not partitioned. To partition this remaining region, we trans- 
late the newly constructed regions downard by the amount A , and continue 
this process until R is completely partitioned. The result of this process 
is illustrated in Fig. 3.3-6 where the entire region -it < 6^ < 0, 

0 < G < Tr/4m is partitioned, for m = 4. Coupling this figure with 
Fig. 3. 3-3 yields the complete .transient behavior of the first order DPLL, 

since the loop responds identically to $ +• 2 tt as it does to 0 . 

o o 

One interesting feature of the digital phase-locked loop is that it 
does not possess an unstable equilibrium point (or, more precisely, se- 
quence), as do both the analog loop and digital loop without VCO harmonics. 
For assume that the VCO is 180° out of phase wi.h the input carrier. 

Then the VCO phase, which is initially tt radians, becomes at k = 8: 


56 




*7 


8 


$ 


T Vi 


-Tr+ 2G 


- C - l)(“2 cospr/8) 
p=C 


and the VCO output is 

w Q = Sc ' 8 tt/S + Q, ) = Sq(2G) = +1 
o c 

Hence the error signal at k = 8 is 

y g = (-2 oosStt/ 3) (+1) * +2 
as opposed to the - 2 originally assumed, and 
P g = -TT+2G + 2G = -tt+4G. 


We see that it is impossible to keep the VCO phase at -tt after eight 
iterations, but instead this phase increases and the DPLL approaches 
the stable equilibrium sequence Cv Fig. 3.2-1. 

A second feature of the DPLL is the variation of the transient time 
as a function of the loop gain, G, for some fixed initial VCO phase. Oi 
c jurse, we expect the transient time to decrease as the loop gain increa- 
ses, as larger loop gains correspondto “faster" loops. This indeed is 
the situation for an analog loop but it is not precisely correct for the digi- 
tal loop, with or without harmonics. 

For the analog loop (where rl.e sum terms are neglected) described by 


we have 


dP 

dt 


- G s in m 


1_ 

G 


0 


dm 
sin P 


A 


P 


O 



(3.3-29) 


Since the left side diverges, we consider the time T for the initial VCO 
phase to reach 0.01 radian- (an arbitrary figure). Then 


58 



* - s\ 


ax 
sia x 


fv 


'0 .01 


and the transient time is inversely g.^portionai to the Joop gain. 

"he digital loop with the VCO harmonics suppressed it described 


by £q. (3.2-1), repeated here for convenience: 


. + 1 = \ ~ (4/")Gsinc? ] 


(3 3-30) 


Here, 0 < ( 4 / rr ) G <2 for stability (Sec. 3.^). The number of iterations 

required for 1$^ [ — 0.01 is obtained as a function of loop gain, G, using 

a simple computer program, and the results are displayed in Fig. 5.3-7 -or 

P - rr/2 and P = -?~/8. The transient t....us are '"■'ir.outec a^ discrete 
o 1 o ‘ 

values of loop gain and so the graph is a set of discrete points. Notice that 
the inverse relationship between gain ar.d i ran s lent time holds only for 
G < 0.S, but not for larger loop gains. From >\pp . 2, we know that for 

0 < (4/tt)G < 1 the VCO phase gees to zero onctonically but that for 

1 < (4 f~) G < 2 it oscillates about zero, its magnitude approach!* ' zero. 
Notice that although the difference equation Eq. (3.3-2) is a discrete version 
uf the differential equation, Eq. (3.3-29), the former exhibits two distinct 
transient behaviors, while the latter exhibits only one. When G is small, 
the difference equation is a good approximation to the differential equation, 
and indeed both equations exhibit similar monotonic transient responses. 
However, if G becomes large*. ( (4/tt)G > 1, to be exact), the approximation 
gets poorer. The difference equation is still stable but now exnioits an osc- 
illatory transient, and the computed results show that the transient incro qsc s 
as the loop gain increases. This is expected as larger loop gain causes a 
greater overshoot of the origin by and hence more iterations are required 
to achieve iP^ [ ^ 0.01. 

For the digital loop, with VCO harmonics induced, the transient times 

are obtained from Figs. 3.3-3 and 3.3-6 by simply following a horizontal 

line at P through the various regions. The .esults aie displayed in Fig. 3.3-8 
o 

for v = n/2 and P = — 3 tt / 8, and we see that the curve is roughly an 
o o 
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Fig. 3.3-7. Transient time as a function of loop gair. for a i*rst order 
D?LL without harmonics . 
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Pig. 3.3-8. Transient time as a function of loop gei:- for a first orcer 
DP LL with harmonics includad. 
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inverse relationship except for certain ir.u«val s or gain where rhc loop 
falls into a "faster* region. 

3.4 .^rsponsc to a Frasue ncv Offset 

If rhr input carrier frequency is deviated by df Hr, a phase error 
develops between the input and YCO, generating a phase error which has 
a nonzero average value. Such a phase error sequence is illustr ate d in 
Fig. 3.4-1. The average value is just the correct amount to deviate the 
VCQ frequency by Hz so that the loop is in equilibrium, tracking the 
deviation with a phase lag error. 

This thesis deals with. the '£>£. demodulation capabilities of a digital 
phase-locked loop as opposed to the tracking capabilities. Therefore, we 
examine die steady state behavior of the D?LL, act the transient response. 
Indeed, a complete solution of Sq. {?. 1-7) is extrei<&iy c;;.~:cuh, even nor 
the case of no raise and an input frequency offset. 

With the DPLZ, in the steady state we can calculate the phase error 
developed as a function of input frequency deviation, and use the result 
to obtain the maximum tradeable offset as a function of the loop gain. 

The input to the D?LL Is 

= -2cos (kTT/2m •> xZuhf/f^) (3.4-1) 

where hf 9 frequency offset in Hz. The VCQ is running at the same fre- 
quency but with a phase error and the VOC output can be written as 

- Sq [(Ic - E)(w/2m + 2*At/f , ) ] 

Let the deviated frequency correspond to 2Jf samples : 

N 9 tTT/(TT/2m + 2wh{/f j 

where the square brackets denote n integer par. . The equation that 

determines E as a function of &f states th- c re.* .C samples the VCO 

phase shall increase by the amount N2nhf/f ; 

s 


Z.**' 


(3.4-3) 


62 



Fig. 3.4-1. Typical phase error sequanoe wr.cn the in,— t frequency is 
dev-atcd from the carrier frequency. 



Pig. 3.4-2. Maxim um tradeable frequency deviation as a function o* 
loop gain for both analog ar.c digital phase-locked loops. 


63 


( 3 . 4 - 4 ) 


S+X-l 


P=£ 


(+1) ! - 2 cos ? (”/2ai + 2nAf/f ) j = N 2^Af/f 


s 


Performing the summation and using the approximation 


N(tt/2m + 2r&f/f ) * n 
s 


we obtain 


(JL . 

2irAf% 

\2m 

TTl 


sin 


-1 


f X2“Lf . I fn _ 2 ~&i 1 I 

; — : sin — i — + — r — : i 


2 L2m f 5 - { , ZA _ 


l 


2G 


> ( 3 . 4 - 5 ) 


J 


The maximum offset trackahle occurs when the argument of the inverse 
sine is unity: 


2G « (N2*r&f/f ) sin(TT/2m+ 2nAf/f )/2 
s s 


(3.4-6) 


Kote that if (n/2m + 2 tt Af/f )/2 « u/2, then Sq. (3.4-6) becomes 

3 


2JrAf «* (4/n)Gf 


s 


(3.4-7) 


which is the corresponding result for a first order analog loop where, 
as observed in Sec 3.1 ( 

G asalog " < 4 /"> 6 f , • P- 4 ' 8 ) 

In Fig. 3.4-2, the maximum trac .e offset, Af, is {dotted against the 
loop gain, G, for both file digital and analog phase-locked loops, for the 
case m * 4, and as anticipated by Bq. (3.4-7), the two curves agree closely. 


3.5 Design Values for the First Order DPLL 

The first order DPLL utilizes a 10-bit A/D co.vvc.~er, and therefore 
the error signal is a 10-bit word. From Table 2 . 5 - 1 , the restriction on 
the phase error signal, , Is 



max 


£ 0.35 


(3 5-1) 
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The res .notion that the first VCQ harmonic fail outsica the D?LL bane- 
width. B,* is given, by Eq. (2.8-4), and repeated here: 

if 

f s /2a - 2 uf > (3.5-2) 

where 

f = sampling frequency 
s 

Af = maximum frequency deviation 

All of the digital phase-locked loops were constructed using TTL and DTL 
logic, having speeds of 1 MHZ and 10 MHz respectively and to allow a 
sufficient computation interval for third order loop calculations, the sam- 
pling frequency is chosen to be 50 KHz. 

The digital filter in the first order DPLL is simply a proportional path 
of gain g^ producing the linearized model of Fig. 3.5-1. If the input 
carrier is deviated in frequency by Af Hz, the input phase is 

<P = 2wAfkT (3.5-3) 

JC s 

and the error signal, e^, appearing in the linearized model of Fig 2.7-1 is 


, 2nAf/f 
4 ' S 


_ _ _ 2rrAf 

*k " n (*/")9 1 «5 vco <3 £ s 


(3.5-4;, 


wtieie G ** ioop gain » g^-Gy^. Imposing the restriction of Eq. (3.5-1) 
yields 

2TrAf 


Gf 


2S 0.35 


(3.5-5) 


The input phase-output voltage transfer function of the linearized 
model is 


IM 

#(z) 


(4/TT) gi ( 1 — Z 


1 • Cl - (4/TT) 0l G yco j 2 


-1 


(3.5-6) 


and taking the difference ( 1-z l ) #(z) as the input ' equency yields the 
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input frequency-output voltage transfer function: 


Xlll. 


(i -*" 1 ) Hz) 




1 - [ 1 - (4/t:}Gj * 


-i 


(3.5-7) 


If (4/~ ) « 1, the bandwidth, 3^, of this transfer characteristic is given 
by (A??. 3) 

» (2/n 2 ) f;G Hz (3.5-8) 


Using Eqs. (3.5-8) and (3.5-5) in Sq. (3.5-2), we obtain the design 
equation 


m ^ f / 11.28 Af 
s ' 


(3.5-S) 


The sampling frequency is chosen to be 50 KHz which allows a 
sufficient computation interval (20 j- sec ) for the DTL and T?L logic 
employed. Then, for Af = 500 Hz, Eq. (2 .5-9) yields 

m ^ 7.4 (3.5-10) 

The VCO algorithm is simplified whan m is a power of two; hence we 
choose m * 4. Eqs. (3.5-5) and (3.5-8) yield 

G ^ 0.216 

B_ *- 2.16 KHz (3.5-11) 

Ai 

The realizable loop gains are of the form 

G » w/5 • 2 N , N = integer (3.5-12) 

• 

and using N * 1 yields G * 0.314, satisfying Eq. (3.5-11). However, 
the maximum allewafcleloop gain from stability considerations is 
it/ 2m * 1^8 * C .394, and *hig proves to be too close to 0.314, so we 
shall use 

G * it/5 • 2^ * 0.15T (3.5-13) 

As the condition ( 4 /tt) G « l is not satisfied, the bandwidth approxima- 
tion of Eq. (3.5-8) is conservative; the actual bandwidth is 
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£_ = 1S00 Hz, 

which still satisfies Eq. (2.5-2) . 

The feet that our maximum frequency deviation is rather small 
( -if = 500 Hz) stems directly from the linitions cf the logic speed. 

Had fester logic been available, then Af cculc have been increased 
by the same factor as the logic speed increases. For example, using 
ECL, the sampling frequency could he increased by a factor of 1 0 and 
hence Af could be ten times larger without changing the lue of m. 
Summarizing, the design values for the first order D?LL are: 

G * 0.15? = it/ 5 • 2 2 

Af = 6C0 Hz 

f * 50 KEz 
s 

m =4 

For an additional illustration of the design procedure, consider a 
voice channel requiring 

f = 3500 Hz 
m 

£ = 3 

Af = 10.5 KHz 

If the A/D converter employs ten bits, we require 

m S f / h. 28 Af = 8.45 * 10 -6 f (3.5-14) 

S' , s 

Using f = 500 KHz allows us to choose m * 4; the required arithmetic 
s 

computations can be performed in one sampling period if ECL logic is 
utilized. The bandwidth restriction is 

^ < f g /2m - 2 Af = 3L KHz 

and loop gain restriction is G ^ 0.3. The choice for the loop gain is 
G*tr/5- j m 0.314 resulting in * 36 XHz. 
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3.«" - rsdv/&re Implementation 

The block diagram of the first order DPLL is shown in Fig. 3.b-l. 
Following the s ample -and -hoi d and analog -to -digital conversion, we 
,.ave a ten-bit binary word, x,^. The bank of exclusive -or gates performs 
the multiplication, of x, and the VCO output, w. = 3q (k“/2m + ? ), 
and the res ult is stored in the output register, whose contents are conver- 
ted back to analog fora and filtered. The output is also scaled and ir.te- 
grated to generate the VCO atgumeet. <k+ l)-/*e * * k+ v is stored 

in tire VOO register and used to multiply the next input, x. . . 

o * 1 * 

To obtain the gain G = ”/5 • 2”, a factor of 1/16 = 1/2* is used for 

the explicit VCO gain, leaving w/ ( 5 /2^) as the implicit VOO gain. This 

implicit is realized by using the third most significant bit of the VCC 

register to feed the exclusive -or gates (multiplier}, as this hit corresponds 
, 2 

to 5/2 volts. Notice that the zowaro Icon contains unity gain; this is 
because the D/A converter uses only tea hits and sacling in the feward 
loop would mean a loss cf bits in the output. Also, the unity gain insures 
that a maximum number of the D/A* s tea hits are used, so that the signai- 
to-guantization noise ratio is a maximum. 

In order to avoid introducing truncation error, the VCO register and 
adder use 14 bits. In this way, all bits are retained after the scaling by 
1/16. The constant ( TT /8)/( /l .25 ) = 1.25./8 has the 14-bit binary 
representation given by 

1.25/8 * 1000010000.0003 

This follows from the fact that the sign bit represents 5 volts. 

To understand the sequence of events in the first order DPLL, consider 
the timing diagram illustrated in Fig. 3.6-2. The computation interval 
(itentical to the sampling period) of 20 u sec is d-v.icd into five addition 
intervals as shown. The last two are for the VCC computation; the remain- 
ing three, labelled ATI, AT 2, ATS are not used «n the rhrst order DPLL — 
they are reserved for the second and third order DPLL computations. Al- 
though Fig. 3.6-1 explicitly shows a three-input adder, only a two-input 
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Hg. 3.6-1. Schematic diii{jroni of the first order DPIJ. 
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Fig. 3.6-2. Timing diagram for the first order DPJ.L. 



accov is employed and we use this two-input adder twice, requiring the 
two VCO addition intervals shown in the timing diagram. 

We start at t = 0 with the input binary word (the A/D converter 
requires 10 nsec for conversion; hence x(t) v/as sampled at t = 10 h sec 
of the previous computation interval). As scon as x. is generated, it is 
passed through the exclusive-or gates; that is, x, is multiplied by 
w, = 3q(kn72m + $,,).♦ At t = 10 is sec, a command is given to store this 
product, y^, in the output register. At t = 13 .usee, is taken from the 
output register and sent to two channels : the output and the feedback 
paths. In the feedback path, (1/16) y is added to the VCO register con- 
tents .denoted by VCO^, and this temporary result is stored in the VCO 
register. During the next VCO computation interval, the contents of the 
VCO register are addea to the cc nstant (Tr/8)/(Tt/1.25) and the result 
stored again in the* VCO register. The VCO register now contains the new 
value of the VCO argument, VCO. , , , where 

VCO kTl = (TT/8)/(TT/i.25) + y k /l6 + VOO k 

This value is used to generate the VCO output by having the third most 
significant bit of the VCO register feed the exclusive-or g^tes. At 
t = 20 nsec, the new input, Xj^ + is available and the computation cycle 
begins again. 

3.7 Output Noise 

In this section we compute the output thermal, quantization, and 
harmonic noise of the digital phase-locked loop. At large input signal- 
to-noise ratios, the DPLL is able to follow the input modulation w..h a 
small phase error; hence the linearized model of £ :c. 2 .8 is used in the 
calculation. This first order linearized model is shov.a in Fig. 3.7-1. The 
quantization noise which arises from tne finite number of voltage levels 
available is present regardless of any thermal noise, ar.d is responsible 
for the maximum attainable output signal-to-noise ratio. As .he quantl- 
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Fig. 3,7-1. Linear! 





zation. noise is present whan the phas^ error is small, the output 
c^XicirviLi zatiou noise is calculated using the linearized model. 

At low input signal -to-noise ratios, the input noise spikes begin 
to appear and the phase error no longer remains small, invalidating 
use of the linearizes model. To gain some insight into the DPLL per- 
formance with input noise spike.., a deterministic model of a noise spike 
is introduced, the DPLL nonlinear difference equation. is solved via 
computer, and the result (VCO phase solution) is examined to see whether 
or not the DPLT f-'Iows the input spikes 

3.7 .1 Output Thermal Noise 

Let the input FM signal be com. led by additive white, Gaussian 
noise having a two-sided power spectral density /2. Then the band-* 
pass filtered signal plus noise is expressed as 

-r - * 

x(t) - UjCtJcos 2Tt f Q t - n 2 (t) sin 2nf o t -r /A cos C2nft+<P(t) j (3.7-1) 

where n^(t) and n 9 (t) are white, Gaussian noises, bandlimited to B^, 
the IF bandwidth.^ end haying two sided power spectral density n. Also, 
i Q is the carrier frequency. When tht input signal-to-noise ratio 
( V N i = T1BJJ./2A 2 ) is large, Eq. (3.7-1) is approximately given by 

x(t) » 2 Acos X 2* ft + <p(t) + n 9 (t)/2A3 (3.7-2) 

That is, the thermal noise introduces a phase contribution 
n(t) = n^(t)/2A, having two-sided power spectral density 

G n * *1/4 A 2 (3.7-3) 

The linearized first order DPLL transfer function is 


H (z) = 


XSsX 

«(z) 


z - 1 

z + G -1 


(3.7-4) 


and the output noise power is 


N 


o 


G (f) | H(e J2nf ^ fs ) | 2 df 
L *m 


(3.7-4) 
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where f is the cut-off frequency cf the cutout low pass filter, aow 
r.i 

f <<: f and so for j f [ ^ f we have 
as s m 


. -• O— £ /r - 

i H(e J ' **) ! * 2^f/Gf, 


(3.7-6) 


y.e.arn: 


2T « (2-r/Gf ) 2 '; 3 / SA 2 

o v - s ' m ' 

2 

The input signal power is S. = 2 A ; the output signal-to- 

<*. 

is 


(3.7-7) 
noise ratio 


K o (2-/Gf ) 2 Tif S /5A 2 

s zr. 


(3.7-8) 


If the modulating signal is c constant, producing a fixed frequency 
deviation if Hz from the earner, the input phase is 


cp = 2Tr£f k /f ^ 


v .i.7-9) 


The linearized model tracks this frequency deviation with a phase error: 


$ v * 9. - 2TrAf/(4/TT)f G 


(3.7-10) 


and the output is a constant: 


y k = 2 trhf /f g r 


(3.7-11) 


Hence the output signal power is S Q = ( 2r r Af/£^ G ) and 


if ■ 3 * 


2 A 

ilf 


m 


(3 7-12) 


as expected. 

If the modulation is sinusoidal .he A put 


usasc ss 


«P, = ]3sin2"i. .. if 

K ill 3 


» 1 * O j 


with 0 ** moduluu*/.. index, f * : adulating frequency, >A, of course. 
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- r = 2 i_ . Than Use linearized loop output is also sinusoidal with 
- 5 ves iy 

<_„ - ° iHCe^^s-,! » 2rgf a /Gf s (3-T-:4) 

2 

The output signal power is S * (y 1 / 2 . end, as expected. 



(2.7-15) 


The results of Eqs. (3.7-12) end (3.7-15) are of course identical to 
reseats for an ordinary discriminator, as the DFLL transfer function is 
identical to a differentiator transfer functioaj that is, above threshold, 
fb* digital phase-leered loop performance is identical to a discriminator. 


^ .i OCfeT-tizatlcn Noise 

Ac ere h*v* a finite somber c£ hits available, we have a finite number 
cz voltage levels representable be these hits. Hence, the sampled input 
signal Is quantized to these levels and an error is com nutted, producing 
a noise component at the output. j3us quantization e component is 
present even with no input thermal noise, and limits the maximum attain- 
able output signal -to-noise ratio. 

A. each sampling instant, when, the A/D converter acts on an input 
voltage level, the binary word obtained represents the voltage level 
rounded off to the nearest available level. Thus the error committed lies 
between -S/2 and + S/2, whe*e S is the quantization step size. Assuming 
that the true level is uniformly distributed over its range, the error is uni- 
formly distributed ever the range (-5/2, S/2), anc has a variance 

a 2 * S 2 / 12 (3.7-16) 

aow, the quartz ation errors from sampi.' .s sample are independent. 
This follows from thu foci that between samples the- input signal varies 
ever many quantization levels. In fact, the sampling h equeccy is below 
the carrier frequency and so, between samples, ihe input varies cv?r at 
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luast one cycle, which means a variation ever all the quantization 
levels. Since the errors are independent, the power spectral density 
of the quantization noise is while [i?j : 


G„ = S 2 / 12 f (3.7-17) 

“q a 


where is the sampling frequency. The octrut quantization noise, 


is ther efo re 
A 1 ® 

3f_ = \ G !H(e 52Sf ^)l 2 df = (2n/G) 2 (f /f ) 3 /lS (3.7-18) 

i*i S 




3JJ Foike Response 

The thermal output noise results obtained above ate valla only et 
large input signal-to-ndse ratios, for which the inpn. noise approxima- 
tion of £q. (3.7-2) and the linearization of the DPLL ere valid. When the 
input signal-to-noise ret'o deteriorates, both of these assumptions ere 
invalid an. the input signal is c or r upte d by noise spile s. A noise spike 
appears when the signal plus noise phasor rotates completely around the 
origin, causing the input phase to increase by 2ft radians. The input fre- 
quency consequently contains a sharp pulse, having the area 2** . Figure 
3.7-2 illustrates the input phase and instantaneous frequency when a spike 


appears. 

The spike duration, T depends upon die IF bandwidth, B_- To 

spuce xr 

obtain this relationship, consider the modulation deviating die carrier to 
the high end of the IF bandwidth. Then s-i the noise components are at a 
lower frequency than die carrier and the tendency f- for the sig r ^u plus 
noise phase angle to rotate clockwise. If at some instant of time all the 
noise power is concentrated at the lowest IF iscuer.cy, then the noise 
phasor is rota&cg at -R^, Hz with respect re- rhe signal phasor, and the 
time required for' one complete i evolution, is 


spike 




(3.7-19) 
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Tzis, if course, represents tile fastest possible spike. 

A sifEsrectiatiOG dL»crittiaator oeaodulates the input phase change 
producing a spike, or, as it is sometimes referred to, a "click", in the 
output. It is tne appoarence of dicks that degrades the output signal- 
end causes the output slgnai-to-noise ratio to fall off sharply. A phase- 
locked loop posse s ses the ability to co nrott an input spike into a doublet 
(the derivative of a spike) by virtue of its periodic phase det e ct or charac- 
teristic. The power spectral density of the doublet Is f* times tire power 
spectral density of the spike and hence the doublet dantritos signifi- 
cantly less power to the baseband signal than does the spike. 

To see hcw a phase-locked loop coovsrts a spike into a doublet, let 
us examine the phase detector characteristic. The characteristic azises 
non the ’Product of input djpfel £ - 2 cos (u Q t tCJ] and VCO square 
wave fundamental [(4/sj sinfurt + $) j and is given by 

« {4/n> sin (* k - (3.T-20) 

Jq. (3.7-20) is plotted in Fig 3.7-3. Wna the phase-locked loop is 
demodulating ahs input, we <ge oparatina in tire vicinity of = 0. 

(Our design equnticn based ecx quai»ti cation noise and haemoaic distor- 
tion requires { j <: 0.35 for a tje-bt A/D converter.) 'When a noise 
spike appears, the input phase rapidly increases by 2^, producing a 
large phase error. Depending on the loop gain, (and the nature of the spike), 
the VCO phase may either Vncrea^S by 2n or remain in the neighborhood 
of zero, m tiie fbnnar case, the DPLL fallows the spike and the output, 
which is proportional to thfc dstivativeof the VCO phase will co nt ai n a .. 
spike, bn the latter cas*,the VCO phase may initially increase (in an 
attempt to follow the spike) but returns to its baitin', value; a typical plot 
of ^ might appear as in Kg. 3.7-4 . Notice th£t contains a spike. 
The output, being proportional to the derivative .. ~e VCO phase, co nt ai n s 
a doublet, as shown in Fig. 3.7-4. Notion, that in order to convert a 
spike into a doublet, the DPIL must '‘slip a cycle" , which is to say that 
the input phase and VCO phase differ by 2n. Tils, of course, has uc 
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effect upon the loop ope r a ti o n as the phase detector characteristic 
is periodic, with period 2iL 

Y/a thse&ra would like to detefioine how the first order DPLL re- 
sponds to input ooise spikes. To this end we introduce a deterministic 
model £or the input phase during a spike: <p fc shall increase from 0 to 
2 tt in T spilca seconds in a sinusoidal manner: 

^spIkaW*^ 1 - 008 ^^^ ° St ^ T spike 

ffinst- ‘ad in Tig. J.T-5a. k addition, we let the input carrier amplitude 
decrease, also in a sinusoidal manner: 

^spike^ * (* + *> -*)/* 3ooa **VT splke (5.7-22) 

for O-St S T^ tt _ .mhare sis foe minimum amplitude. This amplitude 
variation is illnstzatad in Fig. 3.7-55. Bence, when a spike occurs, foe 


*. . -*^v.CM f )OQ«to/8« . » v ♦ 1 (3.7-33) 

A spike is B.OS*: likely to opcar whan foe modulation deviates the 
carrier to one extreme o£ thelP bandwidth, as all foe noise components 
. have a larger for t.uaUex) frequency than thed aviated carrier. We shall 
onnoldior * 

V’. «*•»?»> 

where - : f — widnlatlng frequency 

m 

| w modulation index 

The instantaneous frequency is a maximum Of * Af ) at k * 0; therefore 

v m 

we introduce a negative spike dk *f. In accordance with foe first order 

design values, we choose Af * goo HS and seaside f « 20Q Hz* £ « 3. 

- _ 28 

Then 

Bjp - 2(3+ 1)200 « 1600 He (3.7-24) 
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^tpike ^1 = cos »» * /T spike ^ 



INPUT AMPLITUDE, Ay pa^ tt) = cos 2TCt/ 



Fig. 3.7-5. The deterministic spike model: (a) phase variation; 
(b) amplitude varia ti on . 


"spike * 1 / x60 ° = °- 625 ^ • (3.7-25 

In. terms at the sampling period (T = 20 usee) inis is 

s 

Tspike m * S2 5 /20 U sen = 31 samples (3.7-26) 

the input phase plus the negative spike is illustrated in Fig 3.7-6 for 
. the above parameters* 

Before we examine the DPLL response to the spike, we first examine 
- ■ the DPLL response to the sinusoidal modulation alone so that we have 
\ a basis for comparison. Fig. 3.7«-7 is the VGO phaie, § k> obtained by 
a computer solution c^ foe DFIX di^erence equation: 

\ +1 ~\ ~ ( T */l0)co6(kiT/8 + <P fc ) Sq(krr/8 + 9 k ) (3.7-27) 

with 

9 k = 3 sin 2 ti 200 k/ 50.000 

The solution is generated fay starting at k * -250 with $ ‘ r -« 0, and 
is printed out beginning at t *= 0, so that the transient time is one com- 
. plete cycd§ pftibe n^dylation. The distortion present is attributable to 
two f acts : first, the teleprinter can space only an integral number of 
spaces; second, we ere. looking at the unfntered solution to Eg. (3.7-27) 
which contains the harmonics generated by the square wave VOO. De- 

y 

spite fids *ilstortioa, it is (dear that the DPLL is following the sinusoidal 
modulation. 

The sinusoidal modulation plus noise spike is now introduced to the 
DPLL, and the solution for $ k is found and plotted using the program of 
Fig. 3.7-9. The programming language is FOCAL JlSj} and the program was 
rut on a ?D?-8 minicomputer. Computer variables correspond to equation 
variables according to Table 3.7-1. The program is rather straightfeward: 
lines 2.05 to 2.15 compute the enor signal with or without the spike (the 
proper case being determined by line 2.05), line 2.20 updates the VCQ 
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INPUT PHASE, «p fc 



iig. 5 7-6. The Input sinusoidal phase plus a 32 -sample spike. 
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VCO PHASE, <P k 
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Fig. 3.7-7. First order DPLL respor , to sinusoidal modulation. 
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C- FOCAL* 1969 


01. P5 A T.F 

ffil.’C S tf=3 . 1 A 1 59 J S G=»V5*2*2SS PK=0 
ei.is f K*-25e#i#2/p»;po 2 
iM.se ©01 T 

02.05 S H«3*FSIM<2*W*20e*K/500Sk>>;i C.-i*<T-X»£. 15 

02.10 S EK=<2+F>/2+<2-E>*FC0M2*1«*K/T>/2 

02.11 S KK=-EK*FC0SCX*W/8+H-W*U-FC0S<U*X/i>>> 

02.12 S ETK=EK*FS6'N<FSIX<K*«»/8tFX>>iGOTO 2.2 
02.15 S EK=-R* FCOS < K* W/6+H ) * FSGNE FS I X<K* c /8 +F?1 > 3 
OK. 20 S PK*PX+G*EX 

02. PS t <K 2.9.2.3.2.35 
02.38 F X=0. 1* 60* T 
02.31 T * 

02.35 I <K/10-FIT3< K/10) >2*9* 2.A. 2*9 
P2. AS F X*0* l.FITR<40+<lfl/3>*FK>;T ** *’ 

02. AS T F X*0*1.40;*: " M 

02.50 T **.**. ! 

02.90 C 
* 


Fig. 3.7-8. FOCAL program for solving the first order DPLL equation 
with sinusoidal modulation plus an input noise spike. 
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Program Variable 

Equation Variable 

T 

Spike duration, T spiIce 

E 

Minimum carrier amplitude, e 

W 

Tt 

G 

Loop gain, G ! 

* 

H 

Input phase, 

EK 

Error signal, ; | 

P. 

VOO phase, ' 

: i 


Table 3.7-1. Correspondence between the variables appearing 

in the computer program of Fig. 3.7-8 and the DPIL equation variables. 
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pr'd-'c.. end liner t.30 through 2.90 are the plotting subroutine for 
p’ citing every tenth computed value. 

First, a spike having no amplitude variation ( e= A = 2 ) is intro- 
duced. !The result, shown in Fig. 3.7-9 clearly shows that the DPLL 
follows the input spike, and therefore DPLL output contains a spike, 
if the it amplitude is allowed to drop, the DPLL still follows the spike, 
as illustrated in Fig. 3.7-iO wnere e= 0.5. .Decreasing e further, the 
DPLL suppresses the spike, as evidenced by Fig. 3.7-10 where e= 0.2. 

In general, we expect that as e debases, it becomes easier for the 
DPLL to suppress the input sp.rv., Decause the gem of the DPLL depends 
on trie carrier amplitude. Furthermore, a mop with a small gain cannot re- 
spond quickly to the rapidly cl mcing input phase when a spike is present 


Notice also that when the spike appears, the input anxplitud ' and there- 
fore the loop gain) drops precisely when the phase is changing most rapid- 
ly. Hence, in ch creasing e. we reach a point wh^re the loop becomes 
sluggish enough and -cannot follow the spike. 

If a wider spike ( larger ^ S p^ e ) is introduced we expert the DPLL 
to follow it, as this slower spike will be wit. .in the dynamic capability 
of the loop. For example, let us introduce a 40-semrle spike having 
s= 0.2. Mthough the DPLL div. not follow the corresponding 32-sample 
spike. Fig. 3.7-12 shows that widening the spike allows the ?PLL to 
follow it: in a sense, the DPLL has more time t<^ follow the phase char e. 


even though itf gain may drop. 

In conclusion, tVen, the analysis usi.*g the above noise spike nodei 
shows that we should rot Si*pert any thrashold improvement by the first 
order DPLL mer tne conventional discriminator. It is for the second and 


th*' d order DPLL's where we have more degrees o: tra cm available 
(i.e., additic...^ loop gains to n.djvst) that we expect to obtain threshold 
improvement. 
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VCC PHASE, 



Fig. 3.7-9. First order DPLL response to sinusoidal modulation plus, 
a 32-sample, constant amplitude spike. 
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VCG PHASE, s’r 
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Fig. S - iO. First order DPII. response to sinusoidal rr.oaulation plus 
a 3 2- sample, 0.5 minimum amplitude spike. 
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VCD PHASE, 
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Fig. 3.7-11. First order DPLL response to sinusoidal modulation plus 
a 32-sample, 0.2 minimum amplitude spike. 
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VCO PHASE. £ k 
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Fig. 3.7-12. First order DPLL response to sinusoidal modulation plus 
a 40-sample, 0.2 minimum amplitude spike. 
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5 . ? Karaor-ic Noise 

Sesibes the therrr.al and quantization noise, another noise contribu- 
tion, generated by the VCO square wave harmonics, appears at the output 
of the D?LL. Recall that the sampling frequency is rel ated to the carrier 
via 

{n+l/4a)f s = f o 

wfcera f s “ saa?lin » **<:«**, 

i = carrier frequency 
o 

n, a * integers 

There, ore, the carrier frequency is shifted to the frequency f /4m. The 
VOO square wave contains odd harmonics, which are shifted to the fre- 
quencies 

f /4m, 3f /4m, 5f /4m, .... (2p+l)f /4a, ... 
s s s s 

with p an integer. The phase error signal which is the product cf input 
and VCO signals contains the sum and difference frequencies, iocated.au 

0, f /2m, 2f /2m, 3f /2m, . ,qf /2m. . . . 

s s s s 

with q an integer. Only the term at 0 Hz contains any useful information, 

namely a quantity proportional to the phase difference between input and 

VCO. The design equation developed in Sec. 2.6 guaranteed that the next 

terra, at f /2m Hz, falls outside the DPLL bandwidth. However, since the 
s' 

sampling frequency is equivalent to zero frequency, the term with q « 2m 
is the first to fall into the DPLL bandwidth, and contributes to the output 
noise. This term does not affect the VCO operation (as argued when devel- 
oping the linearized model in Sec. 2.7) but does contribute output. noise 
which we label harmonic distortion. 

To calculate the harmonic distortion, we mus t go back to the VCO and 

find the harmonic terms which generate the distortion. There are two such 
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hsnmonics. the (4m- 1) " and the (4m+ 1) ; the sura of the (4a- i) har- 

monic and the carrier frequency is a distortion term, while the difference of 
the (4a + l) st harmonic and carrier is also a distortion tern. These harmonics 
appear explicitly in the square wave Fourier series as 


Sc(x- = ~ i_ sinx + * • • + sin t -.m -1)xt T ~ ]T sin(4m* I)x + * * 

(3.7-23) 


and v. e have 


-2 cos (xxr/2m + «?.^ ) Sq (krr/2m ■** 9 ^ ) 


— sin (9. -<?,) + 
TT v iC k' 


tr 4m- 1 


sin L(4m- 1 


A 



4m+l 


sin L(4m+ 1)9,, - 9. 

*w A 


+ • 


v v • 4 «*V j 


Hence the distortion terms are 

D “fi£Ti •JaC(*»-l)* !t + «. k 3 

+ tSTI ;s.7-sc; 

The exact amount of distortion observed at the DPLL output (the output 
of the loiv pass filter) depends on the modulation ( 9^ ). We consider the 
case of constant modulation (a frequency deviation from the carrier) and 
sinusoidal modulation. 

If the input signal is deviated from the carrier by Af Hz, then the in- 
put phase tp, is given by 

it 

cp^ * k 2 tt A f/r^ (e.*— •»») 

Using the linearized model of Fig. 3.7-1, we conclude that the first order 
DPLL tracks this deviation with a lag in phase: 
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(3.7-32) 


- <P k - 2wAf / (4/rr)£ s G 


with G, the loop gain. Hence this distortion terms are 

P* ~ j-"-- sinDt2rr4m Af/f g - (4m- l)2nAf/(4/tT)f s G ] 

+ sinCK2iT4m Af/f g - (4m + 1)2^ Af/(4/ir)f s G 3 (3.7-33) 

We observe that when the input is a frequency deviation, Af Hz, the har- 
monic distortion is at the frequency 4m Af Hz. Its amplitude is bounded 
by 


Id I s * 


8m 


TT 16 m 2 - 1 


(3.7-34) 


and its power is bounded by 

»D S20I °9^ ^ TfJtT dB < 3 - 7 - 35 > 

Now the noise in Eq. (3.8-35) will not necessarily be found at the 
low pass filter output since the frequency 4m Af Hz may fall outside 
the filter passband. The filter is set to the maximum modulating frequency, 
f « ( Af) /0, and therefore Eq. (3.7-35) does not apply if 

Af > f m /4m (3.7-36) 

For example, if we are working with f = 200 Hz, ( Af ] = 600 Hz, 

and m « 4, then only frequency deviations less than 200/16 * 12.5 Hz 
contribute harmonic distortion as per Eq. \l.?-35). This is not to say that 
frequency deviations greater than 12.5 Hz contribute no distortion at all; 
they do contribute but not as a result of the 4m ± 1 VCO harmonics. If 
we tabulate the harmonics generated by the multiplier, eventually we find 
one that folds down into the f baseband, and since the VCO harmonic 

IQ 

amplitudes decrease as 1 /n, the noise contribution will certainly be with- 
in the bound of Eq. (3.7-35). 

Next, consider sinusoidal modulation with modulating frequency f 

in 
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and modulation index £ 


<P k - 0 siait 2 nf m /f g . 

According to the linearized loop the VOO phase is sinusoidal : 
\ = Asia(k2Rf m /f s -ra) 


where 


Ae‘ 


J* 


4G/TT 
1 - (1 -4G/1T)2 


-1 


z = e 


J2nf m /f s 


The harmonic distortion terms are 


(3.7-37) 


(3.7-38) 


(3.7-39) 


i 


D = £ 3^3 sin C (4m-l) Asia + a) + 0 sink 2"f m /t s J 


+ ^ 3^3 slnC(4nt«)Asln(lc2>tf n /f s + a) +Pstok2nf n /l s ] (3.7-40) 

Now i m « f g so that in Eq. (3.7-39), z "» 1 and A ** 0 , or « 0 and the 
distortion is approximately given by 


D “ n ie'J- 1 sin (4mp sink 2nf m /f s ) (3.7-41) 

It follows that the harmonic distortion contains odd harmonics of the 
modulating tone. The amplitude of the fundamental frequency (f ) is 
negligible in comparison to the signal amplitude. The (2n+l) s * harmonic 
distortion is 

D 2n« * n ii^Tx 2I 2n+1 (4«^)Stok2n(2n + l)f m /f 8 (3.7-42) 


where J is a Bessel function of order p, and the corresponding noise is 

N « 20101,^^^3 2 J 2 o + .( 4 m « dB < 3 - 7 ‘ 4S > 

2 n+l 

Again, we must be careful and remember that only those harmonic fre- 
quencies passed by the output filter will contribute to the output distortion. 
Eq. (3.7-43) may be bounded by realizing that 4m0 Is large and using 
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(3.7-44) 


asymptotic expression for the Bessel function: 


J U (X) ~ cos (x - nTr/2 - tt/4) 

Hence U n ( x ) 1 5 ^2/wx and 

„ ^ nAt 1 4 8m 2 

Nl) 2rH-l 201og /2 n 16m 2 - 1 /it 2 m £ 


(3.7-45) 


. 3.8 Experimental Results 

The first order DPLL was tested and signal -to-noise curves obtain- 
ed for sinusoidal and constant offset modulation for modulation indices 
of 3 and 10. The parameters used were; 


sampling frequency, f =50 KHz 

s 

loop gain, . G « w / 20 

frequency deviation, Af = 670 Hz 

carrier frequency, f Q = 55.125 KHz (3.8-1) 

The test facility is shewn in Fig 3. S- 1. The carrier amplitude is kept at 
2 volts peak ( input signal power = 3 dB J. 7or constant modulation 
( a frequency deviation from the carrier ) the distortion analyzer is not 
required and the output noise is measured at the filter output. (The Sailan- 
tine meter used does not respond to d.c.; hence its reading dees not include 
the signal power.) 


3.8.1 Generation of the IF Noise 

In order to generate the I? noise, a special technique was employed. 
Ordinarily, one could "rectangular" band-pass filter a white noise source; 
however, it v/us difficult to build such a bandpass filter when 5^= 1300 Hz 
(for 8 = 3 ) or = 1320 Hz (for £ = 10) and f o = 33.125 KHz, requiring 
Q ~ 30. The technique used explicitly generated the r.oise according to its 
bandpass representation 


n(t) « n 1 (t)cos 2nf Q t + » 2 (t)sin2Trf o t 


(3.8-2) 
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Piy. 3.8-1. Tent facility for obtaining signal-to-noise performance curves of the DPLL. 




where n (t) and n. 7 (t) are independent white noise sources bandiimitea 
to 3^,/ 2 Kz and having equal power spectral densities. A block diagram 
illustrating the equipment used appears in Pig 3.8-2. A General Radio 
Gaussian white noise source is filtered using a four-pole filter to ge>~ 
a, (t); n.^t) is obtained similarly. Both n (t) and n (t) arc then mu .died 
by carrier terms and the two products added using a resistive adder. The 
C0° phase shift between the two oscillators is guaranteed by synchroniza- 
tion: the sync out signal of the Wave Teh 132 is precisely 90° out of 
phase with its output signal. When measurements were taken on the DPLL, 
the multiplier outputs were adjusted for equal power output. 

3.8.2 Results for Constant Modulation 

For the case of constant offset modulation the carrier is detuned by 
Af = 6C0 Hz to 53.725 KHz. The output signal is constant and using the 
linearized model, 

y. = 2rrAf/f G = 0.48 volts 

jC S 

which corresponds to the output signal power 
S Q = 10 log = -8.4 dB 

The output quantization noise is calculated from Eq. (3.7-18) and the re- 
sults are: 


j3 = 3, N = -92.4 
< J 

dB 

(3.8-3a) 

0 = 10, N = -107.6 

q 

dB 

(3.8-3b) 


The reason for these extremely low figures lies in the fact that we 

are using a 50 KHz sampling frequency on a signa. bandlimited to 2C0 Hz 

(for £ = 3). As a result of the mammoth oversampling, the quantization 

noise is undetectable. The question one wants to ask, then, is can the 

sampling frequency be reduced ? In fact, reducing f provides more compu- 

s 

tation time. The answer lies in the design equation, Eq. (3.5-9), obtained 
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Gonciation of tin- 








from corn- scoring the VCC harmonics: 
m ^ f^/ii.23 Af 


9) 


Cue must pay a price for reducing f j either Af or m must be reduced. 

3 

hew With i = 50 XHz we are using A; -• SCC Hz ar.d m = 4. Certainly 

we cannot reduce Af significantly and reducing n results in distortion 

contributed by the VCO harmonics. So we use f = 5C l-CHz and conclude 

s 

that the output signal-to-noise rat *. ; will certainly not be limited by the 
quantization noise. 

Figs. 3.S-3 and 3.3-4 are the experimentally obtained results for 

constant modulation with 3 = 3 and 2=10 respectively. The thermal 

2 

noise curve S „/N = 35 S, / *1 £ is superimposed in each case, 
o' o i' n 

The first observation is that both curves level off as anticipated. 

For 3 = 3 the horizontal portion of the curve corresponds to X = -43.5 dB 

G 

while for .S = 10 we have X * -61.5 dB. As mentioned above, this i= not 

o 

due to quantization; it is due to the phase jitter between the carrier oscil- 
lator and the DPLL clock. This jitter was observed or. the sampled carrier, 
whose sampled values are theoretically giver, by X. M Wv X •/ C ^V.iOv.iC 
sequence. The actual sequence obtained is not periodic as the phase jitter 
results in the sampled sequence x *» -2 cos (kn/8 -*•€}, with €, represen- 
tir.g the jitter. 

Using the measured output noise, we may calculate the oscillator 

2 

jitter. Let the phase jitter sequence, e. , have a variance o , indepen- 

•■C ® 

dent of k, and let the phase jitter be independent from sample to sample. 
Then assuming the D/A sequence “looks like" impulses to the output lev/ 
pass filter, the power spectral density of the phase jitter is constant and 
given by 


G = 
e 




s 


with f„ the sampling frequency. The linearizes DPLL has the transfer 


function 


H(f) * 2TTf/f G 

W 
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\w%orc uio approximation is valid lor I f 1 <<: f . The output phase jitter 

m 

noise is 

x j ‘ { " WW s G)\°l/t s )a * H (f)*(y 3 

m 

2 

Using the measured value, X, * -4S.5 dB we find o, = 0.34 radian. 

* 6 

This value is not caused simply by the phase jitter between the carrier 
oscillator and sampling clock, but also represents any input signal dis- 
tortion from a true sinusoid. 

t 

In the laboratory the carrier oscillator was synchronized with the 
sampling frequency in an attempt to reduce the D?LL output noise. The 
synchronization procedure required dividing the sampling frequency by IS 
via a counter, and extracting the seventeenth harmonic of the resulting 
square wave, thus insuring that f Q = { 1 + 1 /IS ) The effect of synchro- 
nization was a reduction in the output noise of 5 dB, which corresponds 
to a phase jitter variance of O.ii radian. As mentioned above, even if it 
were possible to perfectly synchronize the carrier and sampling frequencies, 
any distortion in the sinusoidal nature of the carrier signal introduces an 
effective phase jitter, as the sampled values are in error from their theor- 
etical values. 

The synchronization hardware must be removed when the DP LL ts tested 
with modulation, as it is impossible to both synchronize the carrier and 
modulate it simultaneously. In practice, the carrier frequency Is derived 
using a narrow band, tracking phase-locked loop, and the s ampiing frequen- 
cy derived from it. 

The spectrum of the output jitter was determined experimentally by 
measuring the output power as a function of the 1c .\r pass filter cutoff fre- 
quency. The result was that the output spectrum is parabolic which implies 
e white input spectrum, as the D?Lh differentiates the input phase. 

The second feature of the experimentally obtained curves is the loca- 
tion of threshold. For 0*3, threshold occurs when S,/ tif = 21 dB; for 

* *U 
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p 


u*rG S hOiu Is ov C2s. 


S.S.3 Results for Sinusoidal Mrlulatisn 

The first order DPLL was tested with sinusoidal modulation 

{?. * 2 sin k 2 tt f /f ) with Af*$00 Hz using £ « 3 and £ * 10, 
k ms 

resulting in f = 2C0 Hz and 6G Hz respectively. The output signal 
m 

is sinusoidal with amplitude A and phase a caludated from the DPLL 


transfer function; 


Ae 




£ 


M*. 


-1 


(3.8-4) 


1 + JG 


1 - z 




_ J 2i f *x/^s 


z = e 


For both f « 200 Hz a ad f = CC Hz, we have 

m tu 


e i2ir£ m /J s „ x 


and Sq. (3.8-4) yields 

A » 2*4f/f G * 0.48 volts (3.0-5) 

O 

The output signal power is 

S 0 - l01ogA 2 /2 *** -9.4 dB (3.5-5) 

which is 3 d3 below the signal power for constant offset modulation. 

This is, of course, expected as with constant offset modulation the carrier 
is always deviated by A f, while with sinusoidal modulation the carrier 
deviation is varying from + Af to - Af. 

In the experimental setup the output low pass filter is set to f Hz. 

IC 

Hence the noise measured at the output does not include the third harmonic 
distortion generated by the VQO harmonics. Also, this output filter atten- 
uates the signal power an additional 3 dB so that the experimentally obtain- 
ed signal power is - 9.4 - 3 » - 12.4 dB. Hence we expect the signal-to- 
noise curves to level off 6 dB below the corresponding curves for constant 
offset modulation, as the phase Jitter noise is the same for both typos of 
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Taa experi.xen.tal results are shown in Figs. 3.6-5 and 3.8-6 for 
3=3 and £ = 10 respectively. The curves level off at about 34 dB 
and 43 d3 respectively, which are 6 dB and 7 cB below the corres- 
ponding curves for constant offset modulation, in good agreement with 
the anticipated 6 dB drop. 

The values of* S,/ r tf at which threshold occurs are 18 dB for £ *» 3 
x n* 

and 21 dB for *3 = 10. Threshold for p = 3 agrees closely with that of a 
discriminator [19] as expected from the computer results with simulated 
spikes. Also notice that with sinusoidal modulation, threshold occurs at 
a lower S./ 'If value than with constant offset modulation. This is ex- 

X 9 fll 

pectec as the constant modulation keeps the carrier deviated to its maxi- 
mum value where it is most susceptible to spikes. Sinusoidal modulation 
deviates the carrier to its maximum only momentarily, resulting in fewer 
input spikes at the seme carrier-to-ncise ratio. Hence we have fewer out- 
put spikes with sinusoidal modulation and threshold is improved. 

3.8.4 Results Using a Three-Bit A/D Converter 

The extremely small value for the calculated quantization noise 
prompts us to ask how the quantization step size affects the D ?LL perfor- 
mance. We anticipate that the level portion of the signal -to-nois e curve 
will remain unchanged, as the phase jitter noise due to the incoming sig- 
nal will dominate the quantization noise. If the rest of the curve also 
retrains unchanged, then we can greatly simplify the hardware required 
* v reducing the number of registers, gates, adders, etc, needed, i 

However, decreasing the number of bits deteriorates the phase error 
signal. We are working with a maximum error signal of 0.48 volts; assum- 
ing this signal is uniformly distributed, the error ignal power is - 17.2 dB. 
Using ten bits , the step size is 10/2 10 * 1C mv. and the quantization 
noise is - 50.8 dB. Changing to three bits brings the quantization noise 
up to - 8.8 dB, and the phase error signal is now lost in the quantization 
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Fig. 3.8-5. First order DPLL performance with sinusoidal modulation and (3 
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noise. 

The effect of the number of bits employed was determined exper- 
imentally. The desired number of bits were taken from the A/D output 
and the unused bits tied to the complement of the sign bit. Fig. 3.8-7 
illustrates the DPLL performance using only three bits at the A/D con- 
verter. (this is the smallest number we could use, as two bits would 
provide only levels of ±2.5 v. and 0 v. Our carrier amplitude is 2 v. 

$nd hence all sampled values would be truncated to 0 v. This repre- 
sents a practical, not a theoretical problem.) The step size is now 

3 

S * 10/2 = 1.25 v. and the calculated output quantization noise is in- 

7 2 

creased by the factor (2 ) = 41 dB to -51.4 dB. The quantization 

noise plus the -47 dB jitter noise gives a total output noise of -45.8 dB; 
the measured noise is - 44.8 dB. 

More important than the leveling off noise is the effect of the use of 
fewer bits on the threshold performance. With 3-bit operation, thres hold 
occurs at 24 dB. which is 5 dB above the value for 10-bit operation. 
Furthermore, this value is worse than that for a discriminator, indicating 
that spikes appear at the output when no spikes are present at the input. 
Observe from Fig. 3.8-7 that threshold occurs almost immediately as we 
enter the thermal noise region — aJnost no 45° slope region exists. We 
conclude that the large quantization noise in the phase error signal plus 
the input thermal noise sufficiently perturb the phase error so as to cause 
the DPLL to occas tonally lose lock. The trans lent required to regain lock 
produces a spike in the output, even though no spike was present at the 
input. 
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Fig. 3.JV 7. Fiisl ordor DPLL performance with sinusoidal modulation, p “ 3, and a three-bit A/D converter 




CHAPTER 4 


THE SECOND ORDER DIGITAL PHASE-LOCKED LOOP 

When an integration path is added to the toward loop filter of 
thiu first order DPLL, a second order DPLL results. The additional in- 
tegral path gives us an additional degree of freedom (the integral path 
gain) so that we have more flexibility in designing the second order 
DPLL to meet the design equations of Cb. 3 as well as to produce a 
better spike response than the first ruder DPLL. 

This chapter deals with the properties of a second order DPU., 
including the design procedure for determing die loop gains* the hard- 
ware implementation* and the theoretical and experimental results. 


4.1 Second Order DPLL Structure and Equation 

The block diagram of the second order DPLL appears in Fig. 4.1-1. 

The only difference between this and the first order DPLL is the integral 

path present in the forward loop filter. This integrator is realized in a 

way identical to the VCO integrator design (in fact, the second order 

DPLL was implemented by time-sharing one integrator for both the VCO 

* 

and forward loop filter; the details are presented in Sec. 4.4); the pro- 
portional plus integral filter is shown in Fig. 4.1-2. Representing the 
Z-transforas of e^ and by E(z) and Y(z) respectively, the transfer 
function of this filter is 


mi 

E(z) 


«1 * V 


-1 


1 - z 

and the corresponding difference equation is 

y k * **-1 + ‘'•l + V’lc * 

As with the first order DPLL* we have 
*k * < x k + n k ) Sq(kn/2m+$ k ) 

* -2cos(kTT/2m + cp^) 


(4.1-1) 


(4.1-2) 


(4.1-3a) 

(4»l-3b) 
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ng. Block diagram of the second order digital phase-locked loop 










Vu = § k + G VC0 5 'k > 4 “- Sc > 


w.tere 

e. = error signal 
a 

n, = input If noise 

JC 

9,. = input phase 

/C 

= VCO phase 
K 

G^^q= VCO gain 

Combining Eqs. (4.1-2) and (4.1-3) yields the equc £cr the second 
order DPLL: _ . 

a! t= 29 —'o''*' 

■k*'i k k -1 ; - ;■ 

- (G. + G 9 ) 2 cos (kTr/ 2 m -r 9 . } 3 q (k~/ 2 m $, ) 

+ 2 G, cos C(k - 1 )rr/ 2 r- r?,: . jSq [(k-l)"/ 2 m *r $». . 1 

d. X — 4. X* 1 ! ” 

* (G. t G„ ) m Sq (k~/ 2m -s* $ ) 

a ^ iv X 


-G 2 sq [{k-1 ;~/2m + 9^ _ , } 


X - 1 




where G- = proportional loop gain and ** integral icop gain. As 

. 2 

expected, this is a second c“der, nonlinear diiferen.ee equation in $. . 

- * t • a, 

the VCO phase. 

4,2 Stability Condition on the linearised ICocel 

When the phase difference hatweea the input and VC 3 is small 
( << tt/ 2 ), the phase error is given by 

e k * (4 /it) (-© x - 9 k ) (---2-1) 

and the linearized model of Fig. 4.2-1 me y h; u.,cd. V/e shall u:;e this 
model to find the restriction on the loop yd;..;, G . and C.., tor al.ihiihy. 
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V.5 






The characteristic equation for the linearized model is 


l+ K9i + « 2 -^Kc»:7^r - 0 < 4 - 2 - s > 

i"Z i+z 

which simplifies to 

* 2 + s[ {4/n KGj + G 2 ) - 2 1 -r [ 1 - (4/it)G 1 3=0 (4.2-3) 

with Gj = 9 i^VoO = propor ti onal loop gain 

G 2 “ 9 2* G VGO * 113105X61 loop 9airu 
The D?LL is stable if and only if all of the roots of Eq. (4.2-3) lie in- 
side the unit circle. \z | - 1. The bilinear transformation 

z = (s + l)/(s-l) (4.2-4) 

caps the interior of die unit z-plane circle onto the left-half of the 
s-plane, enabling the application of the Routh-Gurwitz test on. the re- 
sulting polynomial in s. Using Eq. (4.2-4) in Eq. (4.2-3) secures 

(4/w)G 2 s 2 + [•:ht)G 1 s + [4-(4/n)(2G 1 +G 2 )3= 0 (4.2-5) 

for which the Routh-Hurwitz test gives the restrictions for stability: 

(4 /it) G 2 > 0 (4.2-5a) 

(4/n) G. > 0 (4.2-6b) 

• A 

4 - (4/TT)( 2Gj + G 2 ) >0 (4.2-6C) 

The first two conditions are obvious and the third is satisfied when we 
use gains of die form tt/ 5*2 , with N a positive integer. Therefore 
when the phase error Is small, the second order DPLL is stable regard- 
less of how much scaling is Introduced. Of course, using this linear 
mocel, we cannot conclude anything concerning nhe DPLL response to 
arbitrary initial conditions (ini'll integrator value, initial VCO phase), 
as this acquisition problem involves operating initially with large phase 



error, 


4.5 Design of Second Order Loop Gains 

To enable comparison to the first order DPLL. we shall design 
the second order DPLL to demodjlate a maximum frequency offset of 
500 Hz and use m = 4. The two design equations developed in 
Sec* 2*5 and Sec* 2*5 bec ome 

<Vm« S °- SS < 4 - 3 -‘> 

- 1200 > 5, (4.3-2) 

= phase error signal 

= linearized second order DPLL bandwidth 
= sampling frequency = 50 KHz* 

As the second order DPLL trades frequency deviations with zero 
phase error, we consider sinusoidal modulation in connection with the 
condition of Eq. (4.3-1). Let the modulation be at tike frequency f Hz 
with modulation index 0, so that 

<P, * 0sink2rrf /f (4.3-3) 

JC ID S 


where 


f /8 
s' 


\ 

h 


s 


The error signal, e^, is then sinusoidal with amplitude A and phase a 
given by 



(4/n)fe-I) 2 

z* + zt(4/Tr)(G 1 + G 2 ) - 2T + Cl-(4/Tt)G 1 3 


2ee j2tTf n /f 5 


(4.3-4) 

As f « f , we have the approximation z ** 1 in the denominator and 
m s 

I z — 1 |* 2simrf /f ** 2tif ji in the numerator, yielding 
1 m' s m' s 

A “ (2"V* s ) 2 P/G 2 (4.3-5) 
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Hence, condition Eq (4.3-1) is a restriction on the integral loop gain: 

G 2 * 0 ( 2** m /f s ) 2 /0.35 (4.3-6) 

Using £ * 20C Hz, 0 * 3, f * 50 KHz, we have 
XU s 

G 2 ^ O.OC54 (4.3-7) 

We choose 

G 2 * n/5-2 6 “ 0.0098 (4.3-8) 

We cannot satisfy Eq. (4.3-2) as easily because one cannot find 
a simple expression, not even an approximate expression, for the DPLL 
bandwidth, B T , in terms of the loop gains G. and G . . Therefore we 

la 1 A 

calculate the magnitude of the DPLL transfer function 

z(4/tr)(g ; +g 0 ) - (4/TT)g, J 

K(f) - ~2 — : = — - j 

* + z[(4/tt)(G +G <l )-2j+ C1-(4 /tt)G j | 

1 * 1 

(4.S-9) 

as a function of x and determine the bandwidth graphically. In. 

Eq. (4.3-9), G 2 is fixed by Eq. (4.3-8) and G 1 takes on the discrete 
values tt/ 5* 2 N with N an integer. The restriction of Eq. (4.3-2) is 
equivalent to < 4.05 KHz and if we require a factor of two in the 
inequality, we find the possible values of N are N ^ 3. 

* The field of second order DPLL candidates is narrowed by con- 
sidering the response to the simulated noise spike introduced in Sec. 3.7. 
As before we superimpose the spike on sinusoidal modulation, solve 
the DPLL difference equation, Eq. (4.1-1), am. examine the solution to 
determine whether or not the DPLL follows the spike. 

The computer program which achieves this -s shown in Fig. 4.7-1, 
with Table 4.3-1 identifying the variables. This program is essentially 
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f> FOCAL* 1969 
31.05 A N* T* E 

(Ji.13 S tf=3* 1*11595 S Gl=ir/5*2*3«S 62«bV5*£»7;S PK=3i 5 IK=» 

01 . 1 5 ? K=-250* 1*276; DO 2 
3i-20 QUIT 

@2.05 S H«3*F^I»C2*«*2ee*K/500B0>;i CK*<T-K»2.15 
02- IS S EK=<2+E>/2+<£-S>*FC05C2*b*K/7>/2 
02*11 S EX=-EX*FC0SCX*fc/3+H-«*Cl-FC0S<i»*X/T>>; 

62.12 S EK=EK*FSQ«CFSlS<K*U/8+PK»;e0T0 2*2 

02.15 S EK«-2#FC0SCK*»/S+H5*FSG2KFSIN<K*U/3+PK>j 
02.20 S IK=IH>G2*EK;S PX=?X*G1*EK*IK 

32.25 1 <102.9*2.3*2*35 

52.30 F X=0* 1 * 63* T •*." 

02.31 T f 

02.35 I CK/13-FTTRCX/10>>2.9*2.4*2.9 
22.40 F X=e*l*FI?R<4S*<13/3>*PiO;T " ** 

52.45 T "*“*#; F X*0*l*40;7 ** ” 

02.56 T **.'** ! 

22.90 C 
* 


Fig. 4.3-1. FOCAL program to solve the second order DPLL equation 
for the case of sinusoidal modulation plus a spike. 
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» 

1 Computer Variable 

Equation Variable 


i i 

spike duration, T .. 

spuce 

i 

5 

i 

2 

adniicui* ipiLs amplitud' " 

1 

l 

i 

GI 

proportional loop gain, G. 


62 

integral loop gaii , G 2 


H " 

input phase, 

i 

i 

EX 

i 

phase error signal, 

i 

t 

t 

f 

K 

integrator value 

j 

i 

t 

* 

PX 

VCO phase, 

« 

1 

j 


Table 4.3-1. The correspondence between the variables appearing 
in the program of Fig. 4.3-1 and the DP IZ equations. 
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the same as that for the first order DPLL, except for line 2.20 which 
integrates the scaled error signal (G2EK), and includes this value (DC) 
in updating the VCX) phase ( PK). 

The computer results obtained for 3 £ M £ 8 are not very ixa- 
• press! ve. In all cases, when a constant amplitude, 32-sample duration 
spike is introduced, the DPLL loses lode and does not immediately re- 
gain lock. Tig. 4.3-2 is foe VCO phase for the case N =4 and spike 
duration 0.54 msec. It is dear that the DPLL has lost lock in response 
to the spike, and even if lock is finally regained, foe output noise gen- 
erated during the transient period severely deteriorates the DPLL output, 
making it worse than the first order DPLL. Hence it is impossible to 
find a second order DPLL which meets the phase error requirement and 
responds satisfactorily to input noise spikes. 

To continue the search for a second order DPLL requires changing 
foe integral loop gain and therefore changing foe maximum phase error 
signal. The phase error condition was based upon having equal quanti- 
sation noise and harmonic distortion at the phase detector. We have 
already seen that the quantization noise (end hence the error signal har- 
monic distortion) is not significant compared to the phase jitter, and so 
we shall decrease foe integral path by a factor of 2 (to the next availa- 
ble'-gain) and thereby double foe phase error. 

n 

Using Gg * n/ 5 * 2 * t ^ e bandwidth requirement is satisfied for 
N 3. For N * 3, the second DPLL follows a 32-sample constant am- 
plitude spike as shown in Fig. 4.S-3 and therefore we expect this second 
order EPLL to have an identical threshold to the first order DPLL. 

Setting N *4 and introducing a 32-sample, constant amplitude 
spike, we obtain the response of Fig. 4.3-4, in which the DPLL does not 
follow the Input spike. Widening the input spike to 40-samples yields 
the response of Fig. 4.3-5 where the VCO phase initially follows the spike 
but recovers and suppresses the spike. Further widening of the input 
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a 32-sample, constant amplitude spike with N * 3. 
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spike to 4S-samples results in the DPLL reproducing the spike, as 
illustrated in Fig 4.3-6. If the minimum spike sampltude is dropped to 
i volt, the DPLL r.o longer fellows the 49-sample spike; Fig. 4.3-7 
illustrates this. Keeping the minimum amplitude at 1 volt, the spike 
duration must be increased to 70-sampies before the DPLL tracks the 
spike. 

For N = 5, the DPLL follows a 46-sample, constant amplitude spike , 
and follows a 68-sample, 1 volt minimum amplitude spike. 

For N= 6, we find that the DPLL tracks a 40-sample, constant 
amplitude spike. If the minimum spike amplitude is reduced to 1 volt, 
and the spike duration increased to 60 samples, then the DPLL loses lock 
as a result of the input spike. " ‘The response is shown in Fig. cj.4-8, 
where we see it takes two cycles of modulation for the DPLL to regain 
lock. 

Thus the choice for the proportional path gain is between = 4 
and N = 5. In order to eliminate one value we consider the spike model 
superimposed on constant modulation, which results in a 600 Hz fre- 
quency deviation from the carrier. The computer program is modified so 
that H= <+> k = 2nk 600/50000. 

For N = 4 the DPLL responses to various constant amplitude spikes 
are illustrated in Fig. 4.3-S. In (a) and (b) the spike duration is 40 ana 
43 samples respectively and in both cases the DPLL does not follow the 
spike, a 44-sample spike is introduced and the response; (c), shows 
that the DPLL now follows this spike. When the spike is again widened 
to 46 samples, the VCO phase develops an error of 4tt radians (Fig. >.3-Sd) 
’rsfied of the expected 2 tt radians. This occurs because of the DPLL 1 s 
i.r . . .>• : . ‘'•ansient response, which is observed in Fig. 4.3-9c as an 
i. w ji the steady state input phase. In Fig. 4.3-9d, this under- 

shoot becomes sufficient to bring the phase error beyond the unstable e- 
quilibrfum point -jr 3rr, and so steady state is reached at 

ic ” \ “ 4TT * ** eace t* 1 ® DPLL follows the 46-sample spike, and generates 
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Fig. 4.3-6. Second order DPLL response to sinusoidal modulation pluj 
a 49-sample, constant amplitude spike. 
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Fig .4 .3-7. Second order DPLL response to sinusoidal modulation plus 
a 49 -sample, one volt minimum amplitude spike. 
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Fiej. 4.3-8. Second ord'-r DPI.L response to sinsuoidal modulation plus a 60-samplo, one volt minimum 
amplitude spike. 
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Fig. 4.3-9a,b. Second order DFLL response to constant offset modulation 
plus a spike of duration (a) 40 samples, (b) .3 samples. 
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Fig. 4.3-9c,d. Second order DPLL response to constant offset modulation 
plus a spike of duration (c) 44 samples, (d) 46 samples. 
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an output spike having an area of 4u It is possible for the undershoot 
to generate a utt phase error, as illustrated in Fig. 4.3-9e, where the 
input spike lasts 60 samples. Hem the DPLL output contains a 6n-area 
spike in responst tc the input 2u-area spike. 

Next consider N = 5; the spike responses appear in Fig. 4.3-10. 

In (a), we see that the DPLL deos not follow a 32-sample spike. This 
remains the case until the spike duration.. is 38-samples, at which point 
the DPLL follows the spike, shown in (b). Increasing the spike duration 
to 39 samples, the DPLL develops a 4n radian error, again because of the 
undexdamped response, as shown in (c). Note however that this DPLL with 
N = 5 is more undexdamped than the DPLL with N = 4, as here, a 39- 
sample spike generates a 4n-area spike, whereas previously (N * 4) a 
Zrt -area spike was generated. This underdamped behavior causes trouble 
when a 45-sample spike Is introduced; Fig. 4.3-lOd shows that the VCO 
slips several cycles in an attempt to relock to die input, and hence gen- 
erates a large output spike. Xotice that in Fig. 4.3-9d the VCO has not 
yet regained lode and the output spike generated will have an area of at 
least 8 tt. 

The fact that the DPLL having N = 5 is more underdemped than that 
with N = 4 is further borne out by examining the linearized DPLL transfer 
function. The two curves are shown in Fig. 4.3-11 where we immediately 
notice the larger resonant peek for N e 5, Of course, when the phase er- 
ror is large, the linearized model does not apply, but when the DPLL is 
close to its equilibrium value, we can expect a larger overshoot from the 
N - 5 loop than for N * 4. 

Therefore our choice for the second order DPLL gains is 


G x - tt/5*2 4 

(4.3-lOa) 

G 2 - rr/5-2 7 

(4.3-lGb) 
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4.4 Hardware Implementation 

The block diagram of the second order DPLL is shown in Fig. 4.4-1. 

A 

The proportional loop gain G = n /5-2 is obtained using unity gain in 

5 1 

the forward path, a 1/2 scaling for the explicit VOO gain, and an impli- 
cit VCO gain of n/(5/2), obtained by using the second most significant 
bit of the VCO register for the VOO output. This is the maximum implicit 
VCO gain available, as the most significant bit (the sign bit) is always 
1 as a result of adding fire VOO carrier term (ir/2m). The constant offset 
representing xt/ 8 is the binary word 1000100...0. 

7 

& Order to obtain the integral path gain G =tt/ 5*2 , we must add 
* • 

a scaling of 1/2 to the toward loop Integral path. We choose to place 
this scaling before the integration in order to minimize the probability of 
(tvefUmingtte i n teg r a tor. In order to avoid truncation of three bits, fixe 
toward loop i ntegra t o r (adder and register) consists of 13 bits. Hence the 
output adder, in Fig. 4.4-1, is also a 13-bit adder. 

The D/A converter at the output has a 1 0-hit capacity and so a trun- 
cation error is introduced. The output is fed bade to the VOO through a 
5 

1/2 scaling and to avoid truncation, fixe VOO operates using 18 bits. 

As previously mentioned, the toward loop integrator and VOO inte- 
grator are time-shared, which is to say that adders Aj , A^, and 
represent fixe same piece of hardware. Although fixe toward loop integrator 
register contains only 13 bits, it is easily convert e d to an 18-bit word 
by adding five additional hits whose values are the complement of fixe 
sign bit. 

The schematic diagram of the second order DPLL and the associated 
timing diagram appear in Fig. 4.4-2. There are two adders with the addi- 
tional saturation logic and three registers; an output register, a forward 
. loop integral (FLI) register, and the VOO register. The transfer gates pro- 
vide a buffer during the integration process end also provide flexibility 
in the timing. To understand the operation, refer to the schematic and tim- 
ing diagram. The sequence of events Is: 
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Fig. 4.4-1. Block diagram of the second order DPLt 
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Fig. 4.4~2o. Schematic of the second order 13PLL 
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Fig. 4,4' 2b. Second ordor DPLL timing diagram. 



Forward Loop Integration (performed during AT 1): The seeled 
error signal (signal ® ) and FLI register (signal (5) ) are added 
and the result is stored in the FU register. 

Output Addition (AT 1): The error signal and its integral are 
added and stored in the output register. The D/A conversion also takes 
place at this time. 

VCQ Integration (AT 4 and AT 5): First, the scaled output (signal 
(D ) is added to the VCO register (signal @ ) and the result stored in 
the VOO register. This occurs during AT 4. Next, the VCO offset (signal 
(§) ) is added to the VOO register (signal (§) ) and the result is stored 
in the VOO register. 

The second bit of the VOO register is then used to multiply (exclusive-or 

and the sequence begins 

again. 

4.5 Output Noise 

As with the first order DPLL, the second c.-der DPLL output contains 
quantization, thermal, and harmonic noise, hi addition, a truncation noise 
is introduced as the D/A converter uses only ten of the thirteen output 
signal bits. 

4.5.1 Quantization Noise 

The input signal is quantized using B hi ts and a quantization step 
size S, introducing white quantization noise l aving power spectral den- 
sity G w (f) given by 

« 2 

6„ (f) - S' 5 T a /12 

where T is the sampling period. The linearized second order DPLL 
8 

transfer function is 


gate) tee next input (x^ + ^), generating e^ + 
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V' z \ 

* 4 v M / 


[ ^ g i + ^ >J 2 )2 ~f g i ] Cg " 1 3 


m 


z ” z ( !T 9 i G vCO + " 9 2 G VCO ” 2 ^ + ^ ” tt S 1 G VCO^ 

(4.5-1) 


where g^,g 9 , -and G^qq identified in Fig. 4.2-1. For 

z = exp (j2rrf/f ) and f <<: f = maximum modulating frequency, we 
s * in 


have z ** 1 as « f g , and Eq. (4.5-1) becomes approximately 


jH(z) | = 


Y(z) j *, 9 _ f/r f 

5(Z)| j2TTf/f s • 


(4.5-2) 


That is to say, the linearized DPLL transfer function is approximately 
that of a differentiator. The output quantization noise is 


N 


f * ■ 

| iHl* G^CfJdf - (2TtS/G vco ) 2 {f a /f s ! 3 /IS (4.5-3) 


-£m 


4.5.2 Thermal Noise 

Since ; second order DPLL behaves approximately as a differen- 
tiator, the output signal-to-noise ratio at high input signal-to-noise 
ratios is given by the results of Sec. 3.8.1: 

S q /N o = 3 3 2 S i /N j for constant modulation (4 .5 -4a) 

S /N * f-0 2 S, /N. for sinusoidal modulation (4.5-4b) 

O O 2 i i 

Again, these results are identical to those for an. ideal differentiating 
discriminator. 

4.5.3 Truncation Noise 

The effect of truncating the number of bits used is essentially the 
same as quantizing; however, in truncating, we quantize a finite number 
of levels into a subset of these levels. The expression fox 'he cuantiza- 
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tior. noise assumes the error is uniformly distributed. This is not 
strictly true for truncation, as the discrete values give rise to discrete 
errors. 

let the D/A converter operate using N bits and have a step size 
S volts. Assume the binary word available has M + N bits, and the 
last M bits are lost in the D/A conversion process. Then the error com- 
mitted in truncation is 0, =S/2^, - 2S/2^, ...» ±(2^- l)S/2^ volts. 
Furthermore, zero error is made in 2^ of the possible numbers; 

nonzero errors ere each made in 2*^~* numbers. The error probability 
density function is thus a set of impulses, the impulse at the origin hav- 
ing area 1/2^, while the impulses at =pS/2^, 1 ^p ^ 2^- 1, having 
area l/2^ +1 . The variance of this distribution is 


2 M -l 


°„ 2 - S (pS/2 M ) 2 /2 M+1 =■ Y >'4.5-6; 


P=-2 M « 


Note that if M is large, this reduces to the result for an error uniformly 
distributed over the interval [ -S» +S J 

Assuming the truncation error is endependent from sample-to-sample 
and that the D/A staircase looks like a sequence of impulses to the out- 
put low pass filter, the power spectral density of the truncation noise is 
constant [ 17 ]; 



ff (l-2 _M )(l-2’ M " i ) 


(4.3 _ o) 


with T the sampling period. Hence the truncation noise appearing at 
s 

the output is 


N t = | s 2 (f m /f s )( 1 - 2 ' M )( 1 - 2 ’ M ' 1 ) (--.5-7) 

where f is the output filter cutoff frequency. 
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4. ~ Harmonic Noise 

The harmonic noise generated by the square wave VCO is calcu- 
lated in a sinular manner as for the first order D?LL. *“*«o changes that 
must made are: 

1. The VCO pr. use is now calculated using the second order 
D?LL linear model* 

2. The harmonics generated at the error signal are filtered by 
the proportional plus integral filter in the forward path. 

As developed in Sec 3*8.4, the distortion terns are 

+ n totT sin Ct 4:3+1 ) \ - 'k 2 “- 5 ‘ £! 

For the case of constant modulation, the input is deviated from 
the carrier frequency by A f Hz and 

= k2TrAf/f s (4.5-S) 

In the second order DFLL linearized model, the VCO follows this input 
with zero phase error, $. = , and the distortion terms become* 

cV iC 


D » — ~*r "*" sink2Tr4m A-.; /f 

rr l6cr - 1 ' s 


(4.5-10 


As with the first order DPLL, the harmonic distortion occurs at the fre- 
quency 4m A { Hz. The proportional plus integral filter has the transfer 
function 


T(f) 


S 1 + s 2 




(4.5- 


If the harmonic frequency (4m Af ) falls within the output low pass filter 
cutoff (£ m ), then the filter transfer function is approximately 


143 


b^ 



2C * ) U*e ,2r4ia Af/i s 

by Virtue of the fact that 4a Af/f < 

' s 


g 9 f,/2TT4m Af (4.5-12) 

f /f « 1. Hence at the output 

31 S 


the haraonlc distortion amplitude is 


’ D 1 *» 


I6m 2 -1 


Vs 

n 2 m if 


(4.5-13) 


and the output povaer is 


4 5 x 

H n « 20 log ■ d3 (4.5-14) 

/S^lSai -1)« aAf 

When the canler is sinusoidally modulated at f He with module- 

xn 

tioa index 0, we have 


q> =0siak2fff It (4.S-15) 

k n' s 


Using the linear model, the VGO phase is sinusoidal with amplitude 
A and phase a ‘where 


Ae 


P 








(4.5-16) 


and since again we have z ** 1, the amplitude and phase are approximately 
A *» 0, a * 0. Hence, as with the first order DPLL. th ^ distortion is 

D » - — sin (4m0 sink 2nf /f ) (4.5-17) 

16m 2 -! X s 

ard it contains odd harmonics of the modulating frequency. The ampli- 

gf 

tude of the (2n+ 1) harmonic at the output is 
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8m 


‘ 3 Sat-l 1 * (f 


16m 2 -1 




(43-18) 


sad the output noise is 

N d - JO-os I D I //T dB 
2n+l 


(43-19) 


As mentioned in Sec. 3*8.4, only those harmonics below the low pass 
filter cutoff frequency will contribute noise according to 2q. (43-19). 


4.5 Sxpreimectal Results 

The second order DPLL was tested using the following parameters: 


sampling frequency, f 
^ s 

propo rt io n al loop gain, 
integral loop gain, G g 
carrier frequency , i Q 
. modulation Index, 0 
carrier amplitude 


50 KBs 
w/5*2 4 
w/5-2 7 
53.125 SHz 
3, 10 
2 volts 


The test facility was identical to the setup used to test the first order 
DPLL, Fig. 3.9-1, including the generation of the IF noise. 


4.6.1 Results for ponstant Modulation 

With the carrier deviated by Af «■ 609 H2, the DPLL output is 
cons ta nt: 

y k * 2u Af/GygQf g * s 132 volts (4.6-1) 

which corresponds to an output signal power 

8 Q - 20 log y^ « 5.7 dB (4.6-2) 

The output quantization noise computed from Eq. (43-3) is 

fi « 3, N q « - 80.4 dB 

0 « IQ, N q « - 953 dB 
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(4.6-3a) 

(4.6-3b) 


The truncation of three bits by the D/A converter contributes, according 
to ^q. (^3*7) 

* J3 = 3, Ny * -66.6 dB (4.6-4a) 

P « 10, Nj * - 71.8 dB (4 .6 -4b) 

"With, no thermal noise added to the carrier, the output noise measured 
is 

£ * 3, N o = - 35.0 dB (4.6-5a) 

.3 = 1C, X q = - 50.3 dB (4.6-SfcJ 

As was die case with the first order DPLL, the output noise which dom- 
inates die quantization end truncation noise, is caused by the phase jitter 
between the carrier oscillator and sampling frequency dock. Notice that 
ftiii output noise has a parabolic .* pectrua, as narrowing the output filter 

from f = 200 to f = 80 Hz resu-ts in a 14.3 dB drop in seise power: 
ai Ql n 

die theoretical drop is (60/200)* * - 15.3 dB. 

The signal-to-uoise curves for .3 * 3 and 3 * 10 appear in Figs. 
4.6-1 and 4.6-2 respectively. Both curves level off at values dose to 
die corresponding first order DPLL curves. The curves then approach the 
thermal noise asymptotic curves and then fell off sharply as threshold 
appears. The input signal-to-noise values at which threshold oepurs are 

P* 3, » 17 dB (4.6-6 a) 

P«10. S l /nf m « 22 dB (4.6-6b) 

This is a substantial improvement over the first order DPLL performance; 
the threshold extension is 4 dB for P *.3 and 5 d3 for p = 10. 

Oomputer simulations and experimental measurements for a second 
order analog phase-locked loop C 14] have produced the results displayed 
in Table 4.6-1, where the digital phase-locked loop results are included 
for comparison. Notice that both analog and digital phase-locked loops 
have comparable performances for ft “ 3. 
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Fig.v4.6-1. Second order DPLL performance with constant offset modulation and 0 
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Fig. 4.B-2. Second order DPLL performance with constant offset modulation, 




Modu lation Index 

Threshold (S f /*tf . dB) 

4 9ft 


| Analog PIL 

Digital PIL 


l Theoretical J Experimental 

Experimental 

3 

| 16.5 ' 

[ * 

16.5 

17 

12 

21.2 

21.6 

V 

22 0 * 1G) 


Table 4.6*1. C om pa ri son of analog and digital phase-locked loop 
performance. 
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4.6.2 Results for Sinusoidal Modulation 

The second order DPLL was tested with sinusoidal modulation 
with Af «* 600 Hz and modulating frequencies of 200 and 60 Hz 0 * 3 
and 0 * 10 respectively). Using the linear model, the output signal am- 
plitude. .A, is 


* 2 + <£ g i + £<V 2} + < 1 -« g i> 


( 4 . 8 - 7 ) 


s«e 


J2*W f s 


and tits output signal power is 

S 0 * 20IogA//2 dB (4.6-8) 

' The calculated values are 

* 0*3. S Q * 3.54 dB (4.S-Sa) 

0“ 10. S 0 « 2.63 dB (4.6-Sb) 


: The calculated Values of q ^ntizatic»n and truncation are identical 
' totbose for constant modulation calculated previously. As the output 
low pass filter is set at fire modulating frequency, we do not include the 
third harmonic disto r tion. Vflth not thermal noise addeo at the input, the 
paeasunrit output noise is -?■ 



' 0* 3. R * -35 J6 

• - o 

dB 

(4.6-10a) 


0-10, ®r o « - 50.2 

dB 

(4.6- 10b) 


This noise is attributable to the phase jitter between the carrier frequency 
and the sampling frequency and is identical to that measured with con- 
stant modulation. 

The experimentally obtained signal-to-noise curves for sinusoidal 
modulation appear in Figs. 4.6-3 and 4.6-4. Threshold occurs at the 
following input signal-to-noise ratios: 
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Fig. 4.fi- 



17 dB 


(4.6-1 la) 
(4.6-1 lb) 


*• V nf m * 

3*10. S 4 /Tl£ * 18 dB 
1 01 

which represent improvements of 1 dB and 3 dB respectively over the 
first order DPLL. 


4.6.3 nf the Number of Bits 

The second order DPLL was modified so that only three of the ten 
bits from thg-A/D converter were used. As with the first order DPLL. this 

serves to increase the quantisation noise by 42 dB. but the total output 

« 

quantisation noise is still negligible in comparison to the phase jitter 
noise. However, the phase error signal becomes lost in the increased 
quantization noise, and we expect a deterioration of threshold. 

Pig. 4.6-5 is the experimental results for the second order DPLL usir.g 
three A/D converter bits, with sinsoidal modulation and 3 = 3. Notice 
that threshold occurs at an input signal -to-noise ratio of 22 dB. which 
is 5 dB above the threshold obtained using a ten-bit A/D converter. 

. This degradation is very dose to that obtained for the first order DPLL. 


4.6.4 Effort: of Truncation Within the DPLL 

As previously mentioned, a truncation of three bits is introduced 
at die Dfk converter, as it can only accept ten bits. We no determine 
the effect on the DPLL performance of using ten-bit arithmetic through- 
out the DPLL. This means that after scaling, we truncate the result to 
ten bits, introducing truncation noise as shown in Fig. 4.6-6. 

. The noise source N_ corresponds to a truncation of three bits. 

T 1 

while N_ represents a five-bit truncation. From Eq. (4.5-6), the power 
t 2 

spectral densities of these noise sources are 

i-10 




5.5 • 10 watts/Hz 
6.4 • 10 * 10 watts/Hz 


(4.6-12) 
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Using the linearized DPLL model, the transfer functions relating 

N m and N_ to the output are respectively 
• 1 ^2 


H.ffl - ~2 4 — i i 

z +z{- s -r* — r-2)*(i--i-r) 
5*2 5-2 5-2 


z -e J2,,f/f s 


4 4 4 

Z( 5/2 + 5 2 2 } 7 5/2 

S 2 (f) * 1 4 4 4 

« + •(- 4 -T +-H- -*) + (!- “4*) 

5*2 5*2 5-2 


j2T7f/f 
z=e • s 


( 4 . P - 13 ) 

Assuming that these noise sources are independent, the output r.olse is 
f 

. .9 .9 

(4.6-14) 


S m 

(!Hl 


i g n + l h 2 I G *C, ) & 

Tj *2 


m 


Now sii. ~e f « f , we have the approximation for { f { < f : 
ms n 


H l( f) » ?¥Ik 
4/(5-2 7 ) 


H 2 (f) “ 2' 


(4.8-15) 


and the output truncation noise is = -35.9 dB. This is precisely 
the amount of phase Jitter noise measured at the output and therefore 
the output noise should increase to -32.9 dB when ten-bit arithmetic 
is used. 

The above analysis does not provide any information about the ten- 
bit arithmetic on threshold. This result was obtained experimentally; the 
signal-to-noise cun e is displayed in Fig. 4.6- 7. First we observe that 
the curve levels off at 34 dB, which is the anticipated 3 dB below the 
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cu-ve where no internal trur nation is performed. Secondly, we notice 
tl'.w* threshold occurs at an input sig r.al -io-noi.se ratio of 18 dB, which 
is 1 dB above the original value. 




RDER DIGITAL PHASE-LOCKED LOOP 


.1 Third Order DPLL Structure aod Equation 


The third order DPLL is obtained by adding a double integral path 
to the forward loop filter of the second order DPLL, resulting in the struc- 
ture of Fig. 5.1-1. Note that the double integral path is realized by the 
addition of a single integrator following foe existing Integrator of the sec- 
ond-order DPLL. The gain is labeled g^/g^, making the double integral 
path gain (g 3 /g 2 ) 

Denoting the new integrator output by b^* and the previous integra- 
tor output by a^, we have foe equa tion s 

\‘ a k-i +a 2 8 k 3-t-n 

which combine to yield the difference equation of foe proportional plus 
integral plus double integral digital filter: 

“ l *k- 1 ‘ y k-2 + <°1 + 8 * + «S >e k 

- !*«!♦«,)%.! 

+ 8 I e k-2 (51_2i 

In addition, we have for the VOO 

Vi"*k tG vao y k < s - l - s) 


and the phase error 

Sq(kn/2m + ® fc ) (5.1-4] 

Combining Eqs. (5.1-2)* (5.1-3)* and (5.1-4) yields foe difference equation 
of the third order digital phase-locked loop: 
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f ig* 5.1-1,. Structure of the third order DPLL. 




♦ (G| ♦ G 2 ♦ 63 )^ 8 q (*»/*» ♦ $ k ) 

♦ (iGj ♦ o 2 )Vi 8,1 C*- 1 )"/** + ft k-l * 

♦ 6 1 \.j 8q C(k-2)n/2» ♦ 3 (5.1-6) 

* 

Gj * g^GyQQ * proportional loop gala 

G 2 - « 2 ’ <3 VoO * a * te * ral lo0p ® ato 

Gj - gj-GyQQ - doable fetegial loop gala 

»-looe(kw/2» + ® k ) + n k (5.1-6) 



« k - IF note®. 

ttlri outer DKL eqaattoB in a third order, imnUn o or dif ference equatloo 
latte VOO phase, $ k > 

&si r^^vr ?r ««>”*«*» ***** 

When tbe VOO to ctoeefr tracMng the input pha.' am phase error 
Isanradnafedy 

. e^ - -Seos(tor/2»+ * k ) 8q(kff/te + ft fc ) 

* 

* (4/n) (S? k - ^ k ) ♦ harmonic tana (5.5-1) 

ffr>q hK»**g the teens produces the linearized nodal of tbetblrd 

oteDKL. nten h Flo* <4-1. We shall find the veetrictioa 00 the loop 



•» + 










G 

(1 » E " 




or, equivalently, . 


0 


3 


4, 


* * ( n G l + w G 2 + n G 3 “ S) * +(3 * 2 W G l" 


w w 2 


)* 




(5.2-2) 


(5.2-3) 


Stability of tbs linearised aodd is guaranteed If all the characteristic 
roots lie tyMft tbs unit dido la the s— plane. Wo ctto ttfl teasioiaB" 
ttoa 

* - Jri (*-*-«> 

Etioh SBpt ftB interior of the writ c fad t In tin s-plane tin loft hdf 

s-plane. The rwwWng polynomial in s. 


K S 0 * * «®S> •* * < 4 » 6 ! - i G S > » 

♦ (8-«|4l I -*io J -f«5 s ) - 0 (S^-S) 

lssah)eated to a Bmtb-BOwttz test which genera te s the four restrictions: 


s G * * • 

*W S »*S > • 


1 (5.2-6a) 
(b. 2-«b) 

(5.2-6c) 


<*J c **l«»X 4 n 6 l-J G J> > n G *<»- 4 H-*K-^ G »> 


(5.2-6d) 

Hoar bon the h w h ww considerations end the VCO algorithm, the gains are 
restricted bo the Com r/5*2 2 , with K a positive integer. Hence all the 
loop gains ere positive and less than unity so mat conditions (a), (b), and 
le) ans irrMifM . c o nd iti on fd) to • 


133 



(5*2-7) 



|e£ I £ 0*35 (5*3-1) 

f # /2»-2Af £ ^ (5*3-2) 

where 2^ is the b andwidth of tire linearised model. 

As wttn the second order DPLL, we consider sinusoidal modulation 

in' imposing the phase error restriction of Eq. (5*3-1). Let the Input phase 

. • ’* 

be 

* k - 0elak2nf m /f 8 (5^-3) 
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where f is the modulating frequency. Then the phase error appearing in 
in 

the linearized model of Fig. 5.2-1 is sinusoidal with amplitude A and phase 
angle or, where 


Ae^* * £ 


1 / z -i’ 3 

tt ' z 1 • 


♦ »*|n< G i + V G s> - S ] * *[» - n 


J 2 n V f . 


Since f « f * we obtain the approximate result 
m % 8 

A® 1 * - P(»"V f S > S 


(5.3-4) 


(5.3-5) 


Combining this result with the design condition of Eq. (5.3-1) yields 


G s * fHanf^/y 3 /0.35 


(5.3-6) 


Using 3 ■ 3, f * 200 Hz, and f *50 KHz, we have 
81 8 


Gj * 1.2-10 


(5.3-7) 


We shall initially select 


G 3 « 1.54 • 10 « ir/5 2 


12 


(5.3-8) 


As was die case with the second order DPLL, we cannot obtain an 
expression for the linearized model bandwidth. B^, in terms of the loop 
gains. Hence we search the G^-G^ plane for those third order DPLL;s 
which satisfy Eq. (5.3-2), which is 

^ < 5.05 KHz (5.3-9) 


The region of the G.-G. plane that must be searched is restricted by the 

4 m 

discrete values of the gains and the stability condition. Having chosen 


G 4 , the stability condition becomes 
3 



(5.3-10) 
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which is approximately given by 

N + M < 11 . (5.3-11) 

( Eq. (5.3-11) is a sufficient condition for Eq.(5.3-10); we note that N * l f 
M * 11 also satisfies Eq. (5.3-10) . ) A search of this restricted region 
eliminates the values M * 1, 2 as they violate Eq. (5.3-9). 

In order to select one third order DPLL from the set of candidates 
satisfying the design equations, we introduce the spike model and observe 
the various DPLL responses. The third order nonlinear DPLL difference 
equation is solved using the FOCAL program appearing in Fig. 5.3-1; Table - 
5.3-1 identifies the computer variables with the equation variables. This 
program is essentially the same as that used in conjunction with the second 
order DPLL, the only difference being line 2.20, in which the error signal, 

EK, is integrated to give AK, which in turn is integrated to give BK. Then 
the VCO phuse, PK, is updated using the proportional signal, G1 -* EK, the 
integral signal, AK, and the double integral signal, BK. The computations 
allow a transient time equal to one cycle of modulation (250 computations); 
also, the Initial conditions (VCO phase, and both i nte g rat or values) are zero. 

file computer results show that all of the third order DPLL candidates 
lose lode la response to a constant amplitude. 32-sample spike. While we 
have shown that for small phase errors these loops are stable, an input spike 
increases the phase error and the behavior predicted by the linear model is 
no longer valii. It turns out that with Gg « n/5 *2 , the third order DPLL 

cannot regain lock after a large disturbance. 

In an attempt to find a stable third order DPLL, the effect of the double 

13 

integration was reduced by changing Gg to n/5 * 2 . The result again was 

that any combination of Gj and Gg yielded a DPLL that lost lock in response 
to an input spike. 

14 

Further reduction in Gg to n/5 • 2 produces several DPLL* s which 
are stable in response to a constant amplitude, 32-sample duration input 
spike. These candidates were first subjected to constant amplitude, var- 
iable duration spikes to compare their performance with the second order 
DPLL. Only one, with N ** 3, M ■ 8, showed Improvement in that the input 
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C- FOCAL# 1969 
01-10 A L#M#N*T*E 

01-80 5 W»3- 141595 S G1«W/5*8*N5S G2=W/5*2*M5S G3«fc/S*2«L 
01-30 S PK-05S AK"05S BK«0 
01.40 F K«-850» 1 • 6001 DO 8 

01- 50 QUIT 

02- 05 S H«3*FSINC8*V*800*K/50000>5I CK*<T-K»8-15 
08-10 S EK»<2+E>/2+<8-E>*FC0S<2*tf*K/T>/8 

02-11 S EK»EK*FSIN<K*W/8+H-W*t l-FCOS<W*K/T>) > 

08-18 S EK*EK*FSGN< FCOS<K* W/8+PK> ) 5 GOTO 2-8 

08-15 S EK=2# F SI N C K* W/8+H > * FSGNC FCOS<K* B/8+PK»3 

08-80 S AK=AK+68*EK5 S BK=BK+CK*G3/G2J S PK«PK+G1*EK+BK+AK 

08-25 1 (K> 2-9# 2- 3# 8- 35 

08-30 T !!!5F X»0*1#605T '*•** 

08-31 T # 

08-35 X <K/10-FITRCK/10>>2.9#8-4#2.9 
08-40 F X»0> 1 , FI TRC 40+ < 10 /3>*PK> 5 T " •* 

08-45 T •*•**# #1F X«0#1#405T " » 

02-50 T I 
08-90 C 
* 


Fig. 5.3-1. FOCAL program to solve the third Older DP LL equation 
for sinusoidal modulation plus a spike. 
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Program Variable 

Equation Variable 

PK 

VOO phase, 

AK 

Integrator value, a^ 

BK- 

Second integrator value, 

EK 

Error signal, 

G1 

Proportional loop gain, 

G2 

Integra! loop gain, G 2 

G3 

Double integral loop gain, G. 


Table 5.3-1. The correspondence between the variables appearing 
in the program of Fig. 5.3-1 and the third order DPLL equation. 
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spike duration could be increased to 67 samples before the DPLL began 
following the spike. (The second order DPLL began following spikes having 
a 49-sample duration.) Fig. 5.3-2 is the VOO phase in response to a 68- 
sample duration, constant amplitude spike, and it is clear that the DPLL 
follows the spike. If the spike amplitude is permitted to drop to 1 volt, 
the DPLL does not follow this spike as displayed in Fig. 5 3-3. 

When the spike duration is increased to 70 samples, the VOO phase 
responds in an altogether different fashion. Fig. 5.3-4 is the VOO phase 
for a 0.5 volt and 0.2 volt minimum spike amplitude. Notice that now the 
DPLL neither follows the spike nor supresses it; instead, positive spikes 
are generated. In (a), the DPLL relocks to the sinusoida. modulation with 
a + 2 tt phase error and in the piocess generates a + 2 tt - area spike. Note 
that this spike lasts 110 samples (the input spike lasts 70 samples). In 
(b), the DPLL relocks with a 4tt phase error, generating a 4tt - area spike 
lasting 150 samples. What we observe here is that the DPLL loses lock as 
a result of both the phase error being increased and the carrier amplitude 
being decreased by the input spike. In the process of regaining lock, positive 
spikes are generated, and these spikes have a larger duration and area than 
the input spike. 

An attempt to further decrease brings no improvement. While the 
momentary loss of lock phenomenon is not present, the effect of the double 
integral path is so small so that the third order DPLL ne haves essentially as 
the second order loop. For example, with L = 15, M = 7, - 4 3 input 

spike duration can be widened to 47 samples before the . -i.ows it; 
the second order DPLL wit u identical proportional loop gain and integral 
loop gain follows a 49-sample spike. Other gain combinations yield poorer 
performance than the second order DPLL. 

Therefore, we conclude that although a third order DPLL has the abil- 
ity to suppress more input spikes than a sec ad order DPLL, the third order 
DPLL also has the tendency to lose lock in response to wide input spikes. 
While the DPLL regains lock, in doing so it generates spikes whose area 
may be greater than the area of the input spike. 
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pig. 5.3-2. Third order DPLL response to sinusoidal modulation plus a 68-sample, constant amplitude 
spik°. 



VCO PHASE, 4 k 
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F'g. 5.3-3. Third order DPLL resr Mi sinusoidal modulation plus a 
68-sample, 1 volt jninimum amplitude spike. 
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Fly. 5. -8-4o. Third order DPLL response k o »inu30idnl modulation plua a 70-aample, 0.5 volt 
minimum amplitude spike. 
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Fig. 5. 3-4*0. Third oHor DPLL response* to sinusoidal modulation plus a 70-sample, 0.2 volt minimum 
ampll’uclo spl! n. 


5.4 Jantware Implementa t ion 

In older to experimentally verify the predictions obtained via com- 
puter. a third order DPLL having gains 

Gj * rc/5 • 2 3 

G 2 « n/5-2 8 

G 3 * n/5-2 14 

was constructed. 

The Jblock diagram illustrating the gain locations and values appears 
in rig. 5.4-1. where we observe that sce'ing is performed prior to integra- 
tion in order to prevent saturation from occuring in the adders. In order to 
avoid truncation, the registers and associateo adders have increased ca - 
pactty to handle the original ten bits plus foe bits introduced by sading. 
Notice, in fact, that had only ten bits been used throughout, foe DPLL i<ouid 
be Identical to a second order DPLL as foe double integral term is scaled 
(shifted) fay 11 hits, and would not contribute anything to the output. The 
WOO output wh'ch feeds foe exdusive-or aultiolier is foe third most signif- 
■ leant bit introducing the implicit WOO gain of tt/(5/2 2 ). 

Although five adders are shown in Fig. 5.4-1. only two actually appear 
in foe Implemented DPLL- One . s the output adder A^, which employs 21 
bits. The other adder performs foe functions of A . A,. A^. and A^ona tiue- 
shared basis aud employs 26 bits to accommodate the VCO computations. 

The time-shored operation may be understood by referring to Fig. 5.4-2. in 
which a schematic diagram and timing diagram appear. The sequence of e- 
vents is as follows: 

AT j Tlrst integration): The scaled error signal (signal ® ) is 

added to the contents of the forward » '■?. integrator (FLX) register 
signal ® ) and foe result is w in the FLI register, thus 
Integrating the scaled error sigr. i. 

AT 2 (Second Integration): The FLI register cou nts are scaled 
(signal ) and added to the second forward loop Integrator 
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fiQ. 5.4- 2a. Sch«mi'*'r diagram of the third order DPLL, 
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Fig. 5.4 -2b. Timing diagram for tha third ordor DPI.L, 



(FU 2) register contents (signal (3) ) and the result is stored 
in the FU 2 register, thus integrating the integrator value, or 
double integrating the error signc 1 . 

AT g (Sum of both integrations): The FU register contents (signal (§) ) * 

are added to tiie FU 2 register contents (signal (§) ) obtaining 
the sum of integral plus double integral of the error signal. The 
result is not stored anywhere; however, the output register receives 
a command to store the output of the 21 -bit adder, which at this 
time is die sum of the error signal plus the integral plus the double 
Integral of the error signal. Also, the D/A converter is made op* 
erativr , generating the output in analog form. 

AT^(VOO computation): The scaled output (signal (f) ) is added 
to the contents of the VCO register (signal (5) ) and this temp- 
orary sum is stored in the VCO register. 

AT. (VCO computation): The VCO offset (T T /2m = it/8) (signal @ ) 
is added to the contents of the VCO register (signal (§) ) and the 
sum is stored in the VCO register, complettag the VCO phase com- 
putation. The thud bit of the VCO register is used to exdusive-cr 
gate die A/D output, generating the next error signal, and die 

cycle begins again. 

The time-shared operation rakes it a simple matter to change from 
third order to second order operation. All we need do is to move the "load 
output register command pulse " to the second half of AT^. Thai the output 
register contains the sum of only the proportional plus Integral signals, ^nd 
the DPLL operates as the second order system. 

5.5 Output Noise' 

As with the second order DPLL, the third c -d. • DPLL output contains 
quantisation, thermal, truncation, and hann<~ »se. However, the hanr.or.- 
ic noise is now obtained by filtering the multiplier output through a propo - 
tiooal plus integral plus double integral filter. Also, the D/A converter 
truncates eleven bits from the 21-bit output register. 
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5.5.1 Quantization Noise 

As a result of converting the sampled input signal to a B-fait binary 
word, quantization noise is introduced. If the voltages representable by 
the B-bit words are spaced by S volts, and the quantization error is assumed 
to be independent from sample to sample, the quantization noise has a 
power spectral density given by 


G w = S 2 /12 T s (5.5-1) 

where T is the £ ampiing period, 
s 

The linearised DPLL model. Fig. 5.2-1) has the transfer function 
H(z) - Y(s)/ #(z) = 


4 

1? 

[ S.+W** - 1 

t29 l + S 2 )S+9 ll [-■>] 


^« g i +g * +g j >-* 


3-i(2G I+ G 2 ] 

3**<? g i-') 


(5.5-2) 

For s ■ and f << f_, this transfer function approximates to 

(5.5-3) 


s 




whir'i is the transfer function of a differentiator. The output quantization 
notoe la therefore 

f 

N q - £ IhI’g^ df - (inS/G voo ) 2 (f 1I1 /f,) , /l* (5-5-4) 


as was obtained for the second order DPLL, 

Thermal Noise 

Since the third order DPLL behaves a t e . oximately as a differentiator, 
.f/> output signal-to-noise ratio at high inpu. signal -to-notoe ratios is 
Ken by the results of Sec, 3.8-1, which are identical to those of an ideal 
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differentiating discriminator: 





2 



for constant modulation 


for sinusoidal modulation 


(53-5a) 


(53- r /b) 


5.5-3 Truncation Noise 

Tbejresult derived in Sec. 4.5-3 for the output truncation noise 
applies also to the third order Dt*LL : 


Nj - f s 2 <f m /f s )(i - **“)(i-a' M-1 > (s^-s) 

i 

vheie M is the number of bits truncated. 

53-4 Harmonic Noise 

In the third order DPLL, the harmonics generated at Hie phase detec- 
tor are filtered by a proportional plus integral plus double integral digital 
filter having the transfer function T(z): 

T <*> -«, + 9 2 7-47 + 93 77-47-7 < 5 - 5 - 7 > 

I-Z (1-2 ) 

For s * exp (j2uf/f ) with f « f , the transfer function becomes approx- 

8 S 

iioaftety 

T( 2 ) « g 3 / (2nf/f s ) 2 (53-8) 

When constant modulation is introduced, -he carrier frequency is de- 
viated by Af Hz and the third order DPLL trac*. ^ frequency deviation 
with zero phase lag: 

* ^ f / f s (5.5-9) 

The harmonic distortion terms are identical to those n the second order 
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DPLL: 


D * slnk2T;4m i£ /' 3 
and the filtered output amplitude is 


Id i 


4 8m g 3 

tt 16m 2 - ‘ (2n4m &f.f g ) 2 


(5.5-10) 


(5,5-11' 


producing the output noise 

4 8ra g 3 

n D * 20109 7T tt 16m 2 - 1 (2TT4m Af/^ ) 2 dB C>.5-12) 

When the carrier is sinusoidally modulated at f Hz with a modula- 

m 

tioo index ft, the input phase is 

© k = sin k. 2 tt f /f (5.5-13) 

At this frequency, the transfer function from input phase to VOO phase 
in the linearized model is approximately l/_0, and therefore, 

© k ** © k (5.5-14) 

Hence the harmonic distortion terms are identical to those of the second 

order DPLL and contains odd harmonics of the modulating frequency. The 
8t 

(2ml) harmonic ap nears in the output with amplitude 

g 

1 D 2n+1 1 * tt* 16m*- 1 2 *2n+ 1 [2tt (2n+l) f m /f s 3* (5.5-15) 

and contributes n^ise in the amount 


N r 
°2n+l 

This noise component will appear in the output only if the frequency 
(2n+ l)f^ is below the output low pass filter cutofi frequency. 


20 log ( D 2ft+1 ! / / r dB (5.5-13) 
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5.6 Experimental Results 

The third order DPLJ- was tested using the loop gains 

3 

proportional loop gain, G. = tt/5*2 

g 

integral loop gain, G 2 * tt/5 * 2 

14 

double integral loop gain, G_ = tt/5 ♦ ° 

The compute results of Sec. 5.3 show that this set of gains yields a 
stable DPLL which momentarily loses lock in response tc input noise 
spikes. , 

The parameters used in the signal-to-noise measurements are: 

sampling frequency, f - 50 KHz 

s 

carrier frequency, f = 53.125 KHz 

o 

modulating frequency, f = 200 Hz 

s 

modulation index, 0 = 3 

Fig. 5.6-1 la the third order DPLL performance curve for sinusoidal mod- 
ulation. We observe that threshold occurs at an input signai-to-noise 
ratio of 24 dB, which Is 6 dB worse than the first order DPLL. The rea- 
son for this poor performance is found by examining the DPLL output, 
where one observes spikes extending over a full cycle of modulation. 

Such spikes are not simply reproductions of the input spikes but rather 
are the result of the DPLL losing lock. When the input noise is increased 
the output spike durations increase as the DPLL is thrown farther out of 
lode and requires a longer transient to regain lock. 

The above statements are clearly ill-’ -wuted in Fig. 5.6-2 where 
thera appear photographs of the third order DPI J. output when the input 
signal-to-noise ratio is 13 dB. In (a) we obs -v* ordinary noise spikes 
appearing in the sinusoidal, output signal. fb>. the oscilloscope sweep 
speed has been reduced to observe many ^ of the output and we ob- 
serve that the DPLL 1c ;es lock of the sinusoidal modulation and requires 
about six cycles of modulation to relock. Note that in addition, we observe 
tbe ordinary noise spikes occasionally appearing in the output. Finally, in 




(b) 

Fig. 5.6-2. Third order DFLL output with .\.i sinusoidal modulation 
and an input signal -to-noise ratio of 13 dB. In (a), observe input no.<‘ 
spikes being reproduced at the output wHle ii (b, observe the loss of 
lock. 
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(c), a more severe toss of lock is di played, with the DPLL requiring 
about 34 cycles of modulation to regain l^ck. It is precisely this loss 
of lock which contributes significantly to the output noise and is respon- 
sible for the degradation in performance of the third order DPLL. 

The third order DPLL performance with constant offst modulation 
appears in Pig. 5.6-3, where we observe thet threshold occurs at an input 
5 ignal-to-*noise ratio of 17 dB, which is identical: to threshold for the 
second order DPLL. However the third order DPLL occ? icnally loses lock* 
and this accounts for the steeper slope below tn^ ol- then Tor the sec- 
ond order DPLL which retains lock even after an input spike Note for 

example that at S, /df = 15 dB, we have S / N * i9 dB for tht second 
i ' m o / o 

order DPLL and S / N =15 dB for the third order t. ' L. 
o o 

4 computer simulation [14 ] for a third order analog phase-locked 

loop predicts S./'Hf * 15.8 dB at threshold; experimental results yielded 

v m 

Cl4 ] S, /‘Hf = 16.3, which is only a Q r 2 dB improvement over the sec- 
i m 

ond ordei analog phase-locked loop. Notice that the same phenomenon 
occurs with the digital phase-locked loops; the third order DPLL provides 
no threshold improvement over the second order DPLL. 

A pho ograph of an output spike appears in Fig. 5.5-4, for which the 
input signal frequency is de bated 600 Hz above tht. carrier frequency and 
tire input isgnal-to-noise ratio is 14 dB. The spike is negative, as are all 
the spikes, since the inout signal frequency is deviated to its noslfrve ex- 
freme. We may calculate the theoretical area of the outi vt s^xe and com- 
pare it to .he experimentally obtained value. 

Wh m the DPLL follows an ioput noise spik - ?, ihs VCO pnase ir reuses 
by 2tt ra> frens . The DPLL cutput, y^, is related to the VCO phase, q>^* v<a 

Vr\‘ (,,J * 3)I k <5,6 ' 1) 

,'.nd assuming a spike occurs during the Interval k » 1 to 2: * K, the .;ui..nn- 
tion of the output values is 
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Fig. 5.6-3. Third order DPLL performance with constant offset modulation and /3 


Fig. 5.6-4. A thermal noise spike appearing in the third order DPLL 

output when the carrier frequency is deviated by 600 Hz and 

S,/Tlf - 14 dB. 
i in 

I ’5 



Fig. 5.6- •». Illustrating the third order DPLL loss of lock .n response 
to an input tpike. 


188 




-j i‘ k * (40/n >( ® K+i -<P 1# * v )(•'"> ** 80v^rs ~ 6-r) 

As each represents the output voltage over the entire sampling 

period, T * 20 p sec., the area under the output pike i » 

8 

K 

T S y * 1.6 mv-sec (5.6-3) 

8 k-1 K 

Approximating the spike of Fig, 5.6-4 by a triangle, the beneath it 

i 

is approximately 

Area * j(0,8 cm) (3 .6 cm) (0.2 v/cm)(5 msec/cm) * 1 .44 mv-sec. 

(5.6-4) 

in good agreement with the theoretical value. 

The DPLL loss of lock is illustrated in Fig. 5.6-5 in which the . 
oscilloscope sweep is slowed to 50 msec/cm. To put this picture into 
proper perspective .recognize that the two small spikes to the left of the 
large ones are thermal noise spikes similar to that of Fig. 5.6-4. The 
. relative sizes of these spikes clearly demonstrates that the DPLL loss 
of lord: contributes significantly to the output noise. 
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C HAPT ER 6 

CONCLUSIONS 

In this concluding chapter, a summary is presented of what 
has been accomplished by this thesis, and how it compares with 
other phase-locked loop analyses and experiments; in addition, a 
discussion is given concerning what problems have not yet been 
solved. 

The all digital phase-locked loop featured in this dissertation 
was originated by /. Gardnick, D. L. Schilling, and this author, it 
is completely digital in nature; i.e., beyond the A/D converter, all 
signals in the DPLL are binary words. Within the DPLL, there is no 
conversion bade to an analog signal for the purpose of tuning a voltage 
controlled oscillator. In fact, there is no oscillator or counter within 
the DPLL; instead there is an algorithm to determine the value an osc- 
illator would have at the sampling instant. In this way, we generate 
the oscillator output only at the time that its value is needed for other 
computations. 

Another feature of the DPLL presented here is its synchronous, 
real-time operation. The s ampling frequency is constant and all Mie 
required arithmetic and logic operations are performed within one sam- 
pling period, generating an output sequence wheih is converted to ana- 
log form and filtered. An equation (developed in Sec. 2.2) relating the 
sampling frequency to the carrier frequency must be satisfied to guar- 
antee proper DPLL operation. The synchronous operation enables a 
time-shared operation of one DPLL to demodulate several FM signals 
simultaneously. 

To obtain the real-time operation of the DPLL the VCO is design- 
ed so that its output has two possible values, corresponding to a square 
wave oscillator. This eliminates the need for a binary multiplication at 
the phase detector as a multiplication by ±1 is accomplished using ex- 
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ciusive-or gating. While the use of a square wave VCO simplifies 
the DPLL hciawaio, it greatly complicoteo tr»t» analv sis of Li? non- 
linear, DPLL difference equation. A unique feature of this equation 
is that its solution does not converge to zero in response to a carrier 
input but approaches a periodic sequence having an average value of 
zero. This behavior results directly from the absence of any filtering 
immediately after the phase detector multiplier. 

The transient response to an unmodulated carrier could be ob- 

* 

tained only for the first order DPLL and for the case of no input thermal 
* 

noise. The time to gain lock as a function of initial VOO phase and 
loop gain cannot be expressed in closed form but is presented graphi- 
cally. 'When an integrator is added to the forward path (yielding the 
second order DPLL), not even a graphical solution could be developed. 
(Compare to [21 ] where a stability analysis could be performed only 
for a first order system.) As a result, no information about the acqui- 
sition performance of the second and third order DPLL's was obtained. 
No attempt was made to pursue such an investigation. 

The purpose of this work was to design,develop, and analyze a 
digital phase-locked loop for FM demodulation and threshold exten- 
sion. To obtain information about the DPIX performance at low input 
signal -to-noise ratios, a model of an input noise spike was introduced, 
and the DPLL equation was solved using a digital computer. This pro- 
cedure could not predict where threshold occurs, but could only provide 
a comparison between first, second, and third order DPLL' s and a differ- 
entiating discriminator. Despite this limitation, the s pike model was 
successful In finding a second order DPLL which yielded a five dB 
threshold extension beyond that of a first order DPLL. This improvement 
is obtained for constant offset modulation with a modulating index of 10. 
Furthermore, the spike model predicted the momentary loss of lock of 
the third order DPLL in res ponse to an input spike, and it was this loss 
of lock that deteriorated the third order DPLL performance, producing no 
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threshold ext easier* beyonc the second order DPLL for constant 
offset modulation and hitftshoid deterioration for sinusoidal modu- 
lation. This is an extremely important "negative result". In the 
past* investigators have discussed the threshold extension of higher 
order loops. Here* we have shown that threshold extension does not 
occur. 

» 

The experimental results obtained show that the second order 

DPLL provides as much threshold extension as does an analog phase- 

locked loop [l4 ] . Furthermore, it is demonstrated that no additional 
* 

extension is obtained by a third order DPLL* as was previous ly demon- 
strated experimentally with analog phase-locked loops [14 ] . 

The experimental results also shew the effect of the number of 
bits used on the threshold signal -to-noise ratio. It is found that in 
reducing the number of bits in the A/D converter from 10 to 3* threshold 
is deteriorated by 5 dB. In addition, if 10-bit arithmetic is used 
throughout the DPLL ( thus introducing truncation error), threshold 
suffers a deterioration of 1 dB. 

The fact that maximum threshold extension occurs for constant 
offset modulation suggests that further research be conducted in apply- 
ing the DPLL for demodulating M-aiy FSK signals. La this respect, the 
transient behavior of the DPLL must be investigated when the input fre- 
quency deviation aoruptly cnauges. Also* the DPLL acquisition perfor- 
in tne presence ox noise must be examined. 
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APPENDIX 1 


P iRST ORDER D?LL TRANSIENT RESPONSE— 

HARMONICS NOT INCLUDED 

Here we consider the transient response to an unmodulated carrier 
by the first order digital phase-locked loop in which the square wave 
VOO harmonics are suppressed. The equation describing the system is 
Eq. (3.2-1). repeated here: 

§ k+x » \ - (4G/tt) sin $> k (Al-1) 

If 0 < 4G/tt < ^; then any initial except tt (mod *lo 2^ generates a 
sequence whose limit is zero (modulo 2tt). hr general, the sequence has 
infinite length; however, for certain isloated values of (which de- 
pends on the loop gain G ). the sequence has finite length (i.e., the trans- 
ient lasts a finite time). 

Theorem A-l. The sequence generated by Eq. (Al-1) converges to 
zero (modulo 2 *t) if 0 < 4G/w< 2 and n (modulo 2tt ). 

Proof: Observe that, be reuse of the sin ^ term, the sequence is 
unchanged by a 2nn, n an integer, shift; hence we consider -TT< $ 0 < n 
only. 

First restrict the gain G so that 

0 < 4G/rr < 1 

Then the sequence is monotonic decreasing and bounded from below if 

0 < tt, while the sequence is monotonic increasing and bounded 

o 

from above if -tt < $^ < 0. Hence the sequence converges and it is clear 
from Eq. (Al-1) that* the limit is zero, 
ifext, consider 

1 S 4 G/n < 2 

Then the curves (4G/n) sin $ and $ intersect at in the interval 

( -tt, tt), as shown in Fig Al-1. Note that if $ * $ x » then $ x “ 0 { and 
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Fig. Al-1. The determination of 



Fig. Al-2. The determination of i 2< 


A 

*e 
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® 0 for ^ - 2; we reach the steady state after one iteration* 

t A I 

If } <»> o ; < # A , the sequence alternates in sign but decreases in magni- 
tude: 


j <p ( 

k + 1 1 


i® k i 


For when 0 < $ < 4, , we have $ = $ - (4G/rr) sin$ < 0 and 

< $ o since (sin$ o )/$ o < 2/(4 G/n) as (4G/tt) < 2. A sim- 
ilar argument holds when - 4j < < °* Hence the sequence ! 

converges am£from Ed* (Al-1), I I — * 0. Therefore, — ► 0, 

1C 

K ^lio | < rr , there are two regions to consider: 

# I <W O l<# * 

f 2 < l $ o l<n 

where is the VOO phase repuired to generate 4j at the next iteration: 

*1 = # 2 " ( 4 G/ Tr ) stoJ 2 

The value # 9 is shown geometrically in Fig. Al-2. Note that If $ * ± § , 

6 O 6 

then • and ** 0 for k S 2; i.e., the transient requires two 

iterations. 

When 4, < 1$ | < 4 , we have 0< {$ |< 4 , since the curve sin<p 

1 | O o X X 

is between the two straight lines, <P and <P- 4 (see Fig. Al-2). Hence 

V*\ ; 

For i$ I > 4 , we again partition this interval at ±4, where 

O 6 O 


* 2 = * 3 - (4G/TT)sin$ 3 

Now o ** ± 4 generates the sequence * + 4 0 , * * 4 » $, = 0, for 

0 3 X Z Z 1 JC 

k ^ 3 and the transient consists of three iterations. If 4 <|$ |<4 , 

Z O 6 

then I < *2* an< * °* 

This partitioning process is continued inaefinitely yielding the points 

*i • *2 * * * * • V * * ’ ° lven *** 

fj ■ (4 G/tt) sin 
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V > c 

M. — 




If % - $ , then the transient is finite, requiring n iterations, while 
o r 

if tp * , the transient requires an infinite number of iterations, 
o n 
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APPENDIX 2 


. first ORDER DPIL TRANSIENT RESPONSE— 


HARMONICS INCLUDED 


In this appendix we prove the assertions made in Sec 3.3 about the 
first order DPLL transient where the VOO waveform is a square wave. 

Theorem A2-1. If 0 < 2G < n/2m and 


n-l r 

S 2 G cospw/2m + (r - J )n/2m S A < ZJ 2G cos pn/2m + rTr/2m 

p»0 p®0 

then 


w fc * Sq(kn/2m + 


’ +1, 0 ?»k ^ 2m -r 
- 1, 2m -r <k ^ 2m 


(A2-1) 


(A2-2) 


where O^r ^ 2m. 


Proof: We must show that 


0 Skn/2m + $ k <tt for 0^k^2m-r 

TTSkTr/2m + $ k < 2n for 2m-r <k ^ 2m (A2-3) 


These conditions may be simplified by observing that the VOO argument, 
kTT/2m + , is an increasing function of k. This is true because 

* Ge^ ^ -2G > -TT/2m. Hence the conditions of E H . (A2-3) reduce 
to four conditions: 


(i) os$ o 

(ii) (2m-r)rr/2m + $ 2m _ r < tt (0^r^2m) 

(ili) TTS(2m-r+l)n/2m + $ 2nw+1 (0^r^2m) 

(iv) 2mTr/2m + $ 2m < 2 tt 
P roof of these four inequalities is straightforward. 
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(i) From the hypothesis. 


r-1 r- 1 

9 -i — 2Gco^p r */2ir J (> *• I)-’-/2rD ~ 72 (2G cos p^/2r*» + tr/2m) >0 

° p *0 p **0 

(ii' *he VOO operation, 

2m-r-l 

$ 2m-r “ + G 2/ (+1 )( -2G cos pn/ 2m ) 


r 2 m -r-1 

< E 2G cos pTr/2m + rn/2m - E 2GcosptT/2m « 

p=0 p=0 


But 


2 m -r-1 

E cospTr/2m 

p=0 


2 m 

E cos ptr/2m 
p“2m-r 


since 


Also, 


2 m 

E cospn/2m * 0 . 

p=0 


2 m 

E cos ptT/2m 
p=2m-r 


r 

E cos ptT/2m . 

p=0 


Therefore, 

r 2m-r-l 

E cospTr/2ra - E cosprr/2m * 0 

p»0 p=0 

and 

$ 2m-r < r "/ 2m S "• 


(Ui) We have 


A 

© 

2m-r+l 


2m-r 

5 +G S (+l)(-2cospTr/2m) 

p“0 


r-1 2m-r 

^ E 2GcosptT/2m + (r- l)rr/2m - E 2Gcosprr/2m 
p“0 p*0 
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But, as in (li). 


r~l 23i-r 

£ co spn/2m - £ cospn/2m * 0 

p*0 p“0 

and $ ^ (r-l)tr/2m t which is identical to (iii). 

(iv) We have, using Eq. (A2-4), 

2m-r+l 2m- i 

= $ + G £ (+l)(-2cospn/ 9 ra) + G £ (-l)(-2cospn/2m) 

° p=0 p s 2m-r 

% 

r 2m-r+l 2m- 1 

<rTT/2m+£ 2Gcospn/2m- £ 2G cos pn/2m + £ 2Gcosprr/2m 

p=0 p=0 p=2m-r 

2m- 1 2 m- 1 

= rn/2ra + £ 2G cos pn/2m < rn/2m + £ u/2m 

p=2m-r p=2m-r 

= rTr/2m + rTt/2m — 2 n t 
proving (iv) and the theorem. 

Theorem A2-2. With the hypotheses of Th. A2-1, 0 < $ < 

Proof; We have 

2m-r 2m- 1 

+g£ (+l)(-2cospTT/2m) + G £ (-l)(-2cospn/2m) 

p=0 p*2m-r+l 

r-i 

* ^ - 2G - 4G £ cospTt/2m 

p^l 

r-l 

Since £ cospTT/2m > 0, we have immediately $ < 

pal 

Using Eq. (A2-1), 
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r-l r-I 

9, ^ S 2Gcospn/2m + (r-l) Tt/2m - 2G - 4G S cospn/2m 

P=0 p*l 

r-l 

= 2 [ Tr/2m - 2G oos pTr/2m 3 > °« 
p»l 

since 2G < w/2m. 


Theorem A2-3. If 0 £ 2G < n/2n and 

r r-l 

-(r+l)Tt/2m - 2 2Gcospn/2m ^ $ < - rir/2m - 2 2Gcospn/2m (A2-6) 

• P=0 ° p=*0 

i 

* 

and 9 ^ - n, then 

o 

1 -1. 0<kSr 

+ 1. r+l^kS 2a (A2-7) 

where 0 s r ^ 2m- 1. 

Proof: As in Th. A2-1, we need only prove the following: 

(i) 

(ii) rTt/2m + 9 f < 0 
”li) 0 £ (rfl)ir/2m + $ 

(iv) 2mn/2m + $ 2jq < tr . 

The proofs are straightforward: 

(1) is part of*, the hypotheses. 

r-l 

(ii) $ +G 2 (-l)(-2cospn/2m) < -rn/2m, directly from Eq. (A2-6). 

P-0 

r 

(iii) 9 . - 9 + G 2 (-l)(-2cos pn/2m) ^ - (r+l)n/2m, directly from 

* 1 ° p-0 

Eq. (A2-6). 
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r 2m-l 

(iv) « $ +G S(-l)(-2cosp^/2m) + £ (+l)(-2cospir/2m) 

2ra ° p=0 p=i+l 


r-1 r 2 m- 1 

< -rn/2m - 2G £ cos ‘pn/2m + 2G £ cospxr/2m - 2G £ cospn/2m 

p=0 p=0 p=r+l 

fay Eq. (A2-6). Now 


and 


r-1 r 

- 2G £ cospn/2m + 2G £ cos p"/ 2m « 2Gcosrn/2m 

p=0 p=0 

* 2m- 1 r 

- £ cospTr/2m *= £ cos pTT/2m 

p*n-l p=l 


Therefore, 

& < -rTT/2m + 2Gcosrn/2m + 2G 

2m 


£ 

p*i 


cos pn/2m 


= £ C- r /2m + 2Gcospn/2m 1 + 2Gcos m/2rn 
P*1 

Now, since 2G < n/2m, the summation term is negative, so that 

< 2G cost it / 2m ^ n/2a 

proving (iv) and the theorem. 

Theorem A2-4. With the hypotheses of Th. A2-3, $ q < $ 2m < 2G. 


Proof: 


$ + G £ (-l)(-2oospTT/2m) + 


2m 


p*0 


2m- 1 

G £ (+l)(-2cosptT/2m) 

p=rfl 



+ 


2G + 4G £ cos pn/2m 
p-i 
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r 

Since £ cos p"/ 2m > O r we have immediately 
p=l 

ITolng Eq. (A2-6), 

r-1 r 

$ < - rTr/2m - S 2G cos pn/2m + 2G ♦ 4G E oos pTt/2m 

m P*# P*1 

r 

* S C-n/2m + 2GcospT T /2m ] + 2GcosrTT/2m 
< 2Gcos rir/2m ^ 2G. 
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APPENDIX 3 


FIRST ORDER DPLL BANDWIDTH 

In this appendix, we calculate the bandwidth of the linearised 
first order DPLL, having transfer function 

T(s) - * JLai — (A3-1) 

1 - (1 -4G/n)s 

To obtain the frequency characteristec, we set 

z * (A3-2) 

where f is the sampling frequency. The bandwidth, B_ , is the fre- 

S 1* 

quency for which 

|T(e J2WB L/ f S) | » ! T(I) \f/2 (A3-?i 

or, equivalently, 

1 1 - ( 1 - 4G/TT)e~ 12TT ®L/ f s | 2 « 2(4 G/tt ) 2 JA3-4) 

Now, 

| 1 - ( 1 - 4G/ir)e~ 12nBl ^ fs | 2 

- 4(l-»G/n)uin 2 nB T /f + (4G/n) 2 (A3-5) 

la 8 

and if 8G/"« i. then nB^/f g « n/2 and 

4(l-8G/v)sin 2 nB L /f s + (4G/n) 2 » (2" B^ /f # ) 2 + (4G/n ) 2 (*.3-6) 

and we obtain from Eq. (A3 -4) 

" <»/"*) f g G (AS-7) 

The above result has a graphical interpretation. The quantity 
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d - 1 1 — ( 1 -4G/Tr)e~^ e | = |e* 9 -(1 -4G/n) | 

is the distance between z m 1 -4G/n and z * e^ 6 as shown in 
Fig. A3-1. We are looking for 8 such that d */2(4 G/tt). Now if 
4G/n <K 1. then near z * 1, the unit circle and a vertical line passing ' 
through z = 1 are approximately identical, as illustrated in Fig. A3-2. 
Hence to find 8. we may construct an isosceles right triangle as shown, 
and t k e arc subtended by 8 is approximately 4G/n radians. Since 
8* 2 n corresponds to f * f , the bandwidth is 

\ " C f s /* , K«5/ n > - (2/" J )£ a G. 


as before. 

Note that using the vertical line in lieu of the unit circle results 
in a bandwidth which is smaller than the exact value. This is dear 
from Fig. A3-2. 
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z -plane 



Fig. A3-1. Illustrating the quantity d 


|e ,0 -(l-4G/tT)|. 


Izl=l 



Fig. A3-2. The situation in the neighborhood of z * 1 when 
4G/"« 1. 
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II. A New Area of Investigation 

The Phase Lock :-d Loop employs a linear loop filter to process 
the phase detected signal before correcting the frequency of the 
voltage controlled oscillator . 

We feel that the output of the phase detector can be processed 
nonlinearly thereby obtaining a "better" estimate of the phase of 
the incoming signals. Nonlinear processors are in general diffi- 
cult to construct and may even require the use of a digital computer. 
Schilling and Ucci have determined a relatively simple nonlinear 
processor. 

Schilling and Ucci are currently analyzing the response of a 
PLL which uses an adaptive delta modulator as a nonlinear processor. 
The delta modulator can be used to obtain an accurate estimate of the 
phase and is easily integrated, thereby resulting in an efficient low 
volume, low cost device. 

A detailed discussion of this device will b presented in the final 


report 


